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In this paper, we describe the application of two experimental methods on the determi¬ 
nation of the heat release rate during the combustion of wood in stoves. The experimental 
methods are the oxygen consumption method and the carbon dioxide generation method 
and are adopted from the field of fire safety science. After outlining the basic ideas, we 
show the necessary equations and analyze the results from an actual experiment. We 
conclude that both methods appear to be very useful in this application. 

© 2011 Elsevier Ltd. All rights reserved. 


1. Background 

1.1. Wood firing in stoves 

In Finland, approximately 50% of the heating energy for single 
family homes was provided by wood firing in 2007 [1], The 
most common method for utilizing the energy content of 
firewood is to use modern versions of the traditional Finnish 
stove equipped with massive stone or brick walls for heat 
storage. Such stoves are capable of providing a reasonably 
steady heat output even if the heat input strongly varies with 
time, as is inevitably the case during intermittent firing of 
batches of firewood. In fact, a typical mode of operation is to 
fire 1—2 batches of wood each morning; after the burning has 
ceased the chimney damper is closed and the stove requires 
no more attention until next morning. 


In order to study the performance of such stoves, we need 
to know the heat release rate from a batch of firewood burning 
inside the stove. The information can be used for optimization 
of the heat storage capability of the stove and also for the fine- 
tuning of the combustion characteristics of the stove, aiming 
at, e.g., reduction of emissions of harmful compounds. 

1.2. Determination of heat release rate by using 
measurements of mass loss rate 

In principle, the heat release rate could be determined indi¬ 
rectly from mass loss rate measurements. This technique 
requires that the stove be installed on a balance and the 
change of the combined weight of the stove and the batch of 
firewood burning inside the stove be recorded as a function of 
time. Heat release rate is then obtained by multiplying the 
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mass loss rate values by the heating value of the fuel. 
However, this approach creates considerable theoretical and 
practical problems. 

The theoretical problems stem from the fact that the 
burning process is notoriously unsteady and includes the 
partially overlapping phases of heating, drying, devolatiliza¬ 
tion, combustion of volatiles and combustion of char. The 
average heating value of volatile matter evolved from wood 
differs considerably from the average heating value of wood 
char [2,3], Thus, one cannot really estimate heat release rate 
unless additional assumptions are made regarding whether 
the observed weight loss was due to drying, devolatilization or 
char combustion. 

The practical problems stem from two facts. Firstly, stoves 
are generally very heavy when compared with the batch, so 
one has to be able to accurately measure small changes in the 
weight of a heavy object. Secondly, stoves are normally rigidly 
connected to a chimney, which not only is another heavy 
object, but may also transmit loads from other parts of the 
building structure to the stove. Thus, accurate weighing 
becomes quite difficult. 

In view of the difficulties described above, we decided to 
search for other options for the analysis of heat release rate 
during wood firing in stoves. Promising candidates were 
found in the field of fire safety science and are described in 
the next section. 


2. Oxygen consumption method and carbon 

dioxide generation method 

2.1. Fundamental ideas 

The ideas behind the oxygen consumption method and the 
carbon dioxide generation method are illustrated in the 
following. Let C a H b O c be the generic chemical formula for any 
fuel consisting of carbon, hydrogen and oxygen. Assuming 
complete combustion of the fuel to carbon dioxide C0 2 and 
water vapor H 2 0, the combustion reaction can be expressed as 

1 

C„HbO c ± d 0 2 — * aC 0 2 -t- — bH 2 0 (1) 

where 

d = a + ^b-ic. (2) 

Now consider the amount of chemical energy being 
released as heat due to the combustion reaction. Define the 
following heating values: 


Q' = -AH r 

(3) 

AH r 

(4) 

Q =-r 

AH r 

(5) 

Q'" =- - 

a 


where AH r is the heat of reaction for the reaction expressed in 
Equation (1). Thus, Q' is the heating value of the fuel per 1 mol 
of fuel consumed in the reaction; Q" is the heating value of the 


Table 1 - Combustion reactions of selected fuels. The 
heats of reaction are for the case where H 2 0 in the 
products is assumed to be in vapor state. 

Code 

Fuel 

Combustion reaction 

AH r [kj/mol] 

(Rl) 

Carbon 

C 4 0 2 —* C0 2 

-394 

(R2) 

Methane 

CH 4 + 20 2 C0 2 + 2H 2 0 

-800 

(R3) 

Propane 

C 3 H s + 50 2 3C0 2 + 4H 2 0 

-2044 

(R4) 

Benzene 

C 6 H 6 + 7 ^0 2 ^ 6C0 2 +3H 2 0 

-3128 

(R5) 

Methanol 

CH 3 OH + 1 ^0 2 ^ C0 2 +2H 2 0 

-635 

(R6) 

Acetic acid 

CH 3 COOH + 20 2 -> 2C0 2 + 2H 2 0 

-876 

(R7) 

D-glucose 

CgH^Og + 6O2 —* 6CO2 + 6H2O 

-2772 


fuel per 1 mol of oxygen consumed in the reaction; and Q'" is 
the heating value of the fuel per 1 mol of carbon dioxide 
generated in the reaction. The latter two heating values will be 
used in the analysis of the heat release rate by the oxygen 
consumption method and by the carbon dioxide generation 
method, respectively. 

The value of the heat of reaction depends on whether the 
water vapor in the products is assumed to be in liquid form or 
in gaseous form. In this study, the water vapor is assumed to 
be in gaseous form; thus the lower value of the heat of reac¬ 
tion is employed in all calculations, and, consequently, lower 
heating values (LHV) are obtained from Equations (3)—(5). 

Now consider the combustion reactions presented in Table 
1. Here the fuels were chosen to represent a wide variety of 
organic compounds. The heats of reaction were taken from 
[4], except for reactions (R4) and (R6) which are not presented 
in [4], For these two reactions, the heats of reaction were taken 
from [5]. The data listed in [5] are actually heating values per 
unit mass of the fuel; thus, the heating values were first 
multiplied by the molar masses of the respective fuels to 
obtain heating values per 1 mol of the fuel, and the heats of 
reaction were then computed using Equation (3). 

The next step is to apply Equations (3)—(5) to the reactions 
presented in Table 1. The results are compiled in Table 2. 

It can be seen that the value of Q" is approximately 
constant regardless of the fuel. The average value for all seven 
fuels listed in Table 2 is 420 kj/mol, and all values are seen to 
be within ±10% of the average value. There seems to be 
a slight increasing trend as we move from pure carbon to 
hydrocarbons and from hydrocarbons to oxygenated 


Table 2 - Heating values of the fuels listed in Table 1. The 
heating values are for the case where H 2 0 in the products 
is assumed to be in vapor state. 

Code 

Fuel 

Q' 

[kj/mol] 

d 

[] 

Q" 

[kj/mol] 

a 

n 

Q", 

[kj/mol] 

(Rl) 

Carbon 

394 

1 

394 

1 

394 

(R2) 

Methane 

800 

2 

400 

1 

800 

(R3) 

Propane 

2044 

5 

409 

3 

681 

(R4) 

Benzene 

3128 

7j 

417 

6 

521 

(R5) 

Methanol 

635 


423 

1 

635 

(R6) 

Acetic acid 

876 

2 

438 

2 

438 

(R7) 

D-glucose 

2772 

6 

462 

6 

462 

Average: 



420 


562 
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compounds. However, we will neither try to explain nor to 
utilize this trend. 

The fact that the value of Q" is approximately constant 
forms the foundation of the oxygen consumption method: if 
we are able to determine the rate at which oxygen is being 
consumed in the combustion reactions, then, at least for the 
reactions listed in Tables 1 and 2, we can compute the heat 
release rate just by multiplying the oxygen consumption rate 
by the value of Q". In the following sections, we discuss 
whether this observation also holds for wood combustion, 
how to determine the rate at which oxygen is being consumed 
in combustion reactions, and what things must be considered 
in the practical implementation of such experiments. 

As regards Q'", the situation is more tricky. Table 2 shows 
that the values of Q"' vary within a much wider range than the 
values of Q". The highest values of Q'" are for hydrocarbons; as 
we move toward oxygenated organic substances such as 
acetic acid and glucose, the values of Q'" decrease to approx¬ 
imately 450 kj/mol, and for pure carbon we have Q'" = 394 kj/ 
mol. Thus, in principle we could compute the heat release rate 
by multiplying the carbon dioxide generation rate by the value 
of Q"', but in practice the uncertainties are expected to be 
much higher than those associated with the oxygen 
consumption method unless we have at least some idea on 
the composition of the fuel. 

In this work, however, we are here only interested in wood 
firing; thus, the range of variations in fuel composition is ex¬ 
pected to be much narrower than what was displayed in 
Tables 1 and 2. It is therefore quite possible that the carbon 
dioxide generation method will be sufficiently accurate for our 
purposes. In the following sections, we discuss whether this 
really is the case, how to determine the rate at which carbon 
dioxide is being generated in combustion reactions, and what 
things must be considered in the practical implementation of 
such experiments. 

So far, we have only considered the case of complete 
combustion; in particular, all carbon in the fuel is assumed to 
be converted into C0 2 . This assumption appears to be 
reasonable when wood firing in stoves is being studied, since 
high combustion efficiency and low emissions of CO and soot 
are among the main objectives. It may be noted, however, that 
the oxygen consumption method and the carbon dioxide 
generation method were first developed for use in studies of 
accidental fires, which are characterized by low combustion 
efficiency and high emissions of CO and soot. The effect of 
incomplete combustion on the feasibility of the oxygen 
consumption method was therefore analyzed very soon after 
the technique was introduced, and was found to be quite 
insignificant [6]. For the carbon dioxide generation method, 
correction procedures have been proposed to take into 
account the effect of incomplete combustion [7], 

In view of the considerations presented above, we decided 
to launch a study to determine the feasibility of the oxygen 
consumption method and the carbon dioxide generation 
method to analyze wood firing in stoves. 

2.2. Some historical notes 

The invention of the oxygen consumption method is generally 
attributed to Thornton [8], who analyzed data for 121 different 


organic compounds and found that 'the heat of combustion is 
... proportional to the number of oxygen atoms which 
combine with a molecule ... and is found to be equal to 53n 
large calories per gram molecule’. Considering oxygen mole¬ 
cules rather than oxygen atoms and noting that one large 
calorie equals to 4.1855 kj, we obtain 444 kj per mole of oxygen 
consumed in the reaction. Thornton considered higher heats 
of combustion and therefore we should expect slightly lower 
values for lower heats of combustion; this was, in fact, noted 
by Tsuchiya [9] who proposed that when the lower heat of 
combustion is considered, the most appropriate value is 
418.5 kj per mole of oxygen consumed in the reaction. The 
value proposed by Tsuchiya is in good agreement with the 
value obtained here. 

Thornton suggested that his finding might be used for the 
determination of the heat of combustion of fuels of unknown 
composition. We do not know whether anyone ever attemp¬ 
ted this in practice. Thornton did not mention the possibility 
of using his finding for the analysis of the heat release rate 
during burning, which is not surprising at all, since such an 
undertaking requires the possession of an on-line oxygen 
analyzer with great accuracy and negligible drift. Such 
instruments were only developed much later. 

In the field of fire safety science, the application of the 
oxygen consumption method started in the late 1970’s 
[6,9-12], and quickly became an essential part of several 
standardized fire test methods, such as the ISO 9705 Room/ 
Corner test [13] and the EN 13823 Single Burning Item test [14], 
Nowadays the oxygen consumption method is a powerful and 
versatile tool for the determination of heat release rate during 
combustion of organic materials both in fire testing and in 
experimental research of fires. 

In the field of energy engineering, the oxygen consumption 
method has not been used very much. In Finland, however, 
the fuel feeding rate of solid fuel fired boilers has been 
successfully controlled by a technique essentially based on 
the oxygen consumption method [15—17]. In those studies, the 
analysis of the relationship between heat release rate and 
oxygen consumption was carried out in a rather different 
manner with a heavy emphasis on observations of the 
performance of the boiler at different loads, but the data can 
be analyzed to compute the numerical value of Q" for peat 
with the result being 420 kj/mol. In a more recent paper, 
Mazumdar [18] analyzed the relationship between heating 
value and oxygen consumption for different types of coal; 
calculated from his data, the numerical value of Q" is 430 kj/ 
mol. We believe that our study is the first attempt to develop 
a systematic framework for applying oxygen consumption 
method to batch firing of firewood. 

The carbon dioxide generation method is quite analogous 
to the oxygen consumption method. The carbon dioxide 
generation method has also been applied in the field of fire 
safety engineering [5], but not nearly as extensively as the 
oxygen consumption method. As was shown in Table 2, the 
uncertainties associated with the carbon dioxide generation 
method are larger than the uncertainties associated with the 
oxygen consumption method, and this has probably been the 
reason for the greater success of the oxygen consumption 
method. However, our study has rather specific objectives, 
and we felt that the carbon dioxide generation method 
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actually might prove quite sufficient for our purposes; in fact 
a variation of this method was already introduced in an ear¬ 
lier paper [19], Thus, both the oxygen consumption method 
and the carbon dioxide generation method were included 
in our study. 

2.3. Notes on the generality of our approach 

It appears that the average value of 420 kj per mole of oxygen 
consumed in the reaction is applicable not only to the mate¬ 
rials included in Tables 1 and 2 but to a great variety of other 
materials as well. The work of Tsuchiya [9] was already 
mentioned in the previous section. Huggett [6] also analyzed 
a wide range of materials typically forming the fuel in acci¬ 
dental fires in buildings, and recommended a value of 13.1 MJ/ 
kg for the average heating value of conventional organic fuels 
per unit mass of oxygen consumed. This value happens to 
be in perfect agreement with the average value obtained 
here (13.1 MJ/kg can be obtained by dividing 420 kj/mol by the 
molar mass of oxygen, which is 32 kg/kmol). 

The average value proposed by Huggett has been adopted 
to use in international fire test standards, such as [13,14], 
although it is expressed there in relation to the volume of 
oxygen consumed, not to the number of moles nor to the 
mass. Conversion can be obtained by first noting that the 
average molar volume of an ideal gas at a pressure of 101.3 kPa 
and a temperature of 298 K is 0.02445 m 3 /mol. Dividing the 
average value of 420 kj/mol by the average molar volume 
of 0.02445 m 3 /mol yields 17200 kj/m 3 , which is the value 
found in [13,14]. 


3. Objective 

The objective of this study was to find answers to the 
following questions: 

1. Will it be feasible to expand the use of the oxygen 
consumption method and/or carbon dioxide generation 
method into studies of wood burning in stoves? 

2. What, if any, modifications of these methods will be needed 
when compared to their application in fire research? 

3. What uncertainties are expected in the results? 

All these questions are worth some consideration since the 
experimental conditions in studies of wood burning in stoves 
will necessarily be to some extent different from those 
encountered in fire research studies. This point will be elab¬ 
orated in later sections. 


4. Wood as a fuel 

4.1. Chemical composition 

For the purposes of this paper, it was deemed appropriate to 
adopt a simplified description of the chemical composition of 
wood. We assume that wood consists of combustible organic 
material, ash and moisture. The combustible organic material 
is further divided into carbon, hydrogen and oxygen. 


This is not to say that the other elements in wood are 
unimportant; for example, nitrogen content of wood is quite 
important from the point of view of emissions of nitrogen 
oxides. However, the nitrogen content of wood is generally 
less than 0.5% by weight [20,21], and the influence of nitrogen 
oxidation reactions both on the heat release rate from the fuel 
and on the composition of flue gas can safely be neglected for 
the purposes of this study. 

It may also be noted that, in fact, wood does not contain 
ash. Wood contains mineral matter which undergoes a series 
of transformation reactions during the burning of the wood; 
the end product of the transformation reactions is a solid 
residue called ash. The exact nature of these reactions and the 
composition of the residue are not important from the point of 
view of this study, although they can be quite important in 
some other aspects. Thus, we will adopt the commonly used 
term “ash content” to describe the amount of mineral matter 
in wood. The ash content of wood grown in temperate zones is 
generally very small, typically less than 1% by weight [20]. 
Higher ash content is encountered in bark, in some tropical 
woods, and in refuse wood which has been contaminated by 
sand or other extraneous mineral matter, but such fuels are 
outside the scope of our study. For our purposes, it would be 
sufficient to neglect the effect of the ash content of wood 
altogether; however, we will present the subsequent deriva¬ 
tions in a general form in which the ash content is also taken 
into account. 

The moisture content of wood varies considerably; for 
freshly cut wood it is between 30 and 50% by weight, and after 
air drying for approximately one year, it is reduced to between 
18 and 25% by weight [22]. This is the moisture content nor¬ 
mally aimed at in the utilization of firewood in stoves. Note 
that here the moisture content is defined as the mass of 
moisture divided by the total mass of combustible organic 
material, ash and moisture. At this level, the moisture content 
is expected to have an appreciable effect both on heat release 
during combustion and on the composition and quantity of 
flue gas. The effect of moisture will therefore be analyzed in 
this study. It may also be noted that wet bark and other refuse 
fuels having a moisture content in excess of 50% are also 
burned at the furnaces of paper mills and other wood pro¬ 
cessing plants, but such applications are outside the scope of 
this study. 

The proportions of carbon, hydrogen and oxygen atoms in 
wood are not arbitrary but vary within a rather narrow range. 
We will preserve full flexibility in expressing wood composi¬ 
tion by adopting the chemical formula (CH 2m O„) y , and our 
knowledge of the composition of wood suggests that the 
values of both m and n will be slightly less than but reasonably 
close to unity. Tables 3 and 4 show the chemical composition 
of dry wood for a variety of wood species. 

In both tables, the first section shows the mass fractions of 
carbon, hydrogen, oxygen and ash. These form the primary 
information that are normally reported for fuels and were also 
given in the sources used in this study. The mass fractions do 
not always sum up exactly to unity, which is in some cases 
due to round-off errors, and in other cases due to the fact that 
small quantities of sulfur, nitrogen and/or other elements 
have also been reported in the original data. The deviations 
from unity are in all cases insignificant. Throughout this 
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Table 3 — Chemical composition of dry wood for a variety of wood species [22]. 


Wood species Chemical composition of dry wood 



[] 

U>H 

[] 

[] 

W Ash 

[] 

Sum 

[] 

m 

n 

n 

n 

m — n 

n 

m + n 

[] 

Cedar, white 

0.488 

0.064 

0.445 

0.004 

1.000 

0.783 

0.683 

0.100 

1.466 

Cypress 

0.550 

0.065 

0.381 

0.004 

1.000 

0.714 

0.519 

0.194 

1.233 

Fir, Douglas 

0.523 

0.063 

0.405 

0.008 

0.999 

0.723 

0.581 

0.142 

1.304 

Hemlock, western 

0.504 

0.058 

0.414 

0.022 

0.998 

0.690 

0.616 

0.074 

1.307 

Pine, pitch 

0.590 

0.072 

0.327 

0.011 

1.000 

0.731 

0.415 

0.316 

1.147 

Pine, white 

0.526 

0.061 

0.413 

0.001 

1.000 

0.694 

0.589 

0.105 

1.283 

Pine, yellow 

0.526 

0.070 

0.401 

0.003 

1.000 

0.801 

0.571 

0.229 

1.372 

Redwood 

0.535 

0.059 

0.403 

0.002 

0.999 

0.662 

0.565 

0.097 

1.227 

Ash, white 

0.497 

0.069 

0.430 

0.003 

1.000 

0.836 

0.649 

0.187 

1.485 

Beech 

0.516 

0.063 

0.415 

0.007 

1.001 

0.727 

0.602 

0.125 

1.329 

Birch, white 

0.498 

0.065 

0.435 

0.003 

1.000 

0.782 

0.655 

0.128 

1.437 

Elm 

0.504 

0.066 

0.423 

0.007 

1.000 

0.783 

0.631 

0.152 

1.414 

Hickory 

0.497 

0.065 

0.431 

0.007 

1.000 

0.784 

0.651 

0.133 

1.435 

Maple 

0.506 

0.060 

0.417 

0.014 

0.998 

0.713 

0.618 

0.095 

1.331 

Oak, black 

0.488 

0.061 

0.450 

0.002 

1.001 

0.749 

0.692 

0.058 

1.441 

Oak, red 

0.495 

0.066 

0.437 

0.002 

1.001 

0.803 

0.663 

0.140 

1.465 

Oak, white 

0.504 

0.066 

0.427 

0.002 

1.000 

0.784 

0.635 

0.149 

1.419 

Poplar 

0.516 

0.063 

0.415 

0.007 

1.001 

0.727 

0.602 

0.125 

1.329 


[ ] means that a quantity is dimensionless (e.g. a mole fraction). 


paper, the superscript ** is used to indicate quantities that are 
expressed relative to the mass of dry wood. 

The second section shows the values of m and n, which are 
used in this study to express the chemical composition of the 
wood, as was noted earlier, and also the values ofm-n and 
m + n, the importance of which will become apparent in the 
mathematical analyses of the oxygen consumption method 
and the carbon dioxide generation method. 

The values of m and n were computed from the equations 

m = 6 ( 6 ) 

u;J* 

and 


It is worth noting that the variations in the chemical 
composition of different wood species are quite small. In 
particular, it can be observed that all values of m - n seem to 
be within ±0.2 of the average value of 0.1, and all values of 
m + n seem to be within ±0.2 of the average value of 1.3. This 
observation will be utilized in the later analyses. 

For future needs, we adopt a “standard composition” of 
wood, for which uij* = 0.52, tuff = 0.06 and uiq = 0.42. For 
such wood, m = 0.7 and n = 0.6. Here the mass fractions of 
carbon, hydrogen and oxygen sum up to unity even 
without the ash, which is quite acceptable, since the mass 
fraction of ash in wood is typically less than 0.01. As will be 
shown later, the only thing that really matters is that the 
values of m - n and m ± n should be realistic; for the 
“standard composition”, m - n = 0.1 and m ± n = 1.3 as 
expected. 


Table 4 - Chemical composition of dry wood for a variety of wood species [23]. 


Wood species 




Chemical composition of dry wood 



[] 

[] 

u>o 

[] 

W Ash 

N 

Sum 

[] 

m 

n 

n 

n 

m — n 

n 

m + n 

[] 

Birch 

74 

0.498 

0.065 

0.434 

0.003 

1.000 

0.783 

0.654 

0.130 

1.437 

Corsica pine (waste) 

119 

0.473 

0.057 

0.413 

0.056 

0.999 

0.723 

0.655 

0.068 

1.378 

Lodgepole pine (waste) 

124 

0.511 

0.059 

0.382 

0.047 

0.999 

0.693 

0.561 

0.132 

1.253 

Pine sawdust 

136 

0.531 

0.060 

0.393 

0.013 

0.997 

0.678 

0.555 

0.123 

1.233 

Site a spruce (waste) 

181 

0.490 

0.058 

0.430 

0.021 

0.999 

0.710 

0.658 

0.052 

1.368 

White pine 

304 

0.526 

0.061 

0.412 

0.001 

1.000 

0.696 

0.587 

0.108 

1.283 

Oriental spruce 

982 

0.519 

0.061 

0.409 

0.005 

0.994 

0.705 

0.591 

0.114 

1.296 

Pine sawdust 

1123 

0.510 

0.060 

0.428 

0.001 

0.999 

0.706 

0.629 

0.076 

1.335 

Ponderosa pine 

1339 

0.493 

0.060 

0.444 

0.003 

1.000 

0.730 

0.675 

0.055 

1.406 

White fir 

1439 

0.490 

0.060 

0.448 

0.003 

1.001 

0.735 

0.686 

0.049 

1.420 

Pine chips 

1786 

0.513 

0.047 

0.403 

0.030 

0.993 

0.550 

0.589 

-0.039 

1.139 

Spruce 

2220 

0.473 

0.063 

0.461 

0.002 

0.999 

0.799 

0.731 

0.068 

1.530 


[ ] means that a quantity is dimensionless (e.g. a mole fraction). 
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4.2. Heating values 


This section contains an analysis of the heating value of dry 
wood. The effect of the moisture content in wood is postponed 
to a Section 7. The reason is that the effects of moisture extend 
not only to heating value but to flue gas composition as well, 
and the discussion of flue gas composition will only be pre¬ 
sented later. 

Let q DF be the lower heating value of dry wood. The 
subscript DF derives from “Dry Fuel” and was chosen to indi¬ 
cate that we prefer a more general term than dry wood (DW). 
The value of q DF can be determined using standardized labo¬ 
ratory methods, such as, e.g., [24], although this is usually 
limited to research and development activities. In the everyday 
practice of wood firing, the value of q DF is normally not known. 

To assess the feasibility of the oxygen consumption 
method, we need to know the value of Q£ F . This is the amount 
of heat released during the combustion of dry wood divided by 
the number of moles of oxygen consumed in the reaction. 
Similarly, to assess the feasibility of the carbon dioxide 
generation method, we need to know the value of Q£ F . This is 
the amount of heat released during the combustion of dry 
wood divided by the number of moles of carbon dioxide 
generated in the reaction. The values of Q" and Q'" for some 
pure compounds were already presented in Table 2, but the 
fact that wood is a mixture of many compounds means that 
now we wish to start from q DF , not from the heats of reaction. 
Again, it may be reminded that we are considering lower 
heating values, that is, all the reaction products, including 
water vapor, are assumed to exist in gaseous form. 

The number of oxygen atoms consumed during the 
complete combustion of 1 kg dry wood is 


ANo 2 



1 

l + 2( m 


n) 


( 8 ) 


where a negative sign was introduced to indicate that oxygen 
is being consumed in the reaction. Thus, the amount of heat 
released during the combustion of dry wood divided by the 
number of moles of oxygen consumed in the reaction is 


Qd 


F 


^DF 


1 + 2 ( m 


M c qDF 



(9) 


The number of carbon dioxide atoms generated during the 
complete combustion of 1 kg dry wood is 


ANeo 2 = N" 



( 10 ) 


Thus, the amount of heat released during the combustion of 
dry wood divided by the number of moles of carbon dioxide 
generated in the reaction is 


Q& = 


q DF M c qDF 


AN; 


co 2 


l + -(m - n) 


Q£f 


( 11 ) 


The values of q DF , Q£ F and Q£' F for the wood species listed in 
Tables 3 and 4 are presented in Tables 5 and 6. 

It is worth noting that the values of Q£ F seem to be 
approximately constant. For the wood species listed in Tables 5 
and 6, all values of Q£ F are within ±10% of the average value of 


Table 5 - Lower heating values of dry wood for the wood 
species listed in Table 3 [22]. 

Wood species 

Heating values of dry wood 

q** 

[kj/kg] 

Qdf 

[kj/mol] 

Q£'f 

[kj/mol] 

Cedar, white 

18113 

424.22 

445.41 

Cypress 

21495 

427.63 

469.16 

Fir, Douglas 

19631 

420.58 

450.43 

Hemlock, western 

18723 

429.79 

445.78 

Pine, pitch 

24710 

434.05 

502.58 

Pine, white 

19341 

419.53 

441.65 

Pine, yellow 

20781 

425.30 

474.09 

Redwood 

19708 

421.66 

442.05 

Ash, white 

19201 

423.70 

463.32 

Beech 

18974 

414.91 

440.91 

Birch, white 

18668 

423.10 

450.10 

Elm 

19014 

421.11 

453.16 

Hickory 

18712 

423.87 

452.06 

Maple 

18594 

420.62 

440.62 

Oak, black 

17668 

422.49 

434.64 

Oak, red 

18741 

424.75 

454.43 

Oak, white 

19017 

421.14 

452.42 

Poplar 

19344 

423.00 

449.51 


425 kj/mol. The average value of Q£ F for wood appears to be 
quite similar to the average value of Q" for the chemical reac¬ 
tions presented in Tables 1 and 2. This observation increases 
our confidence on the feasibility of the oxygen consumption 
method and we propose that the value Q£ F = 425 kj/mol be 
adopted for subsequent analyses. 

As regards Q^ F , all values are between 411 kj/mol and 
474 kj/mol with the single exception of pitch pine, for which 
Q£' F = 503 kj/mol. Based on what was discussed in Section 2.1, 
we conclude that such a high value indicates hydrocarbon- 
type chemical composition. Table 3 indeed shows an excep¬ 
tionally low oxygen content for pitch pine, which observation 
is, after all, not surprising considering the name of this 
particular wood. Thus, we propose that the value Q(( F = 446 kj/ 
mol be adopted for subsequent analyses. This value is 
approximately in the middle of the range observed for the 


Table 6 - Lower heating values of dry wood for the wood 
species listed in Table 4 [23]. 


Wood species Heating values of dry wood 


q“ Q£f Q£'f 

[kj/kg] [kj/mol] [kj/mol] 


Birch 

74 

18700 

423.20 

450.60 

Corsica pine (waste) 

119 

16930 

415.35 

429.51 

Lodgepole pine (waste) 

124 

17510 

385.72 

411.19 

Pine sawdust 

136 

19450 

414.11 

439.55 

Sitca spruce (waste) 

181 

18250 

435.60 

446.94 

White pine 

304 

19370 

419.19 

441.90 

Oriental spruce 

982 

18760 

410.33 

433.76 

Pine sawdust 

1123 

19030 

431.27 

447.76 

Ponderosa pine 

1339 

18710 

443.28 

455.42 

White fir 

1439 

18650 

445.82 

456.73 

Pine chips 

1786 

19570 

466.99 

457.78 

Spruce 

2220 

17930 

439.89 

454.88 
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other wood species and also satisfies Equation (11) when we 
set Q^p = 425 kj/mol, m = 0.7 and n = 0.6. 

Finally, it should be noted that the data presented in [22] are 
the higher heating values, and the lower heating values used 
in this study were therefore computed using the equation 


<2df — <1df,h — 


Mh,o 

m^7 


w^L — cJdf.h 


9uj]j*L 


( 12 ) 


where q D F,H is the higher heating value and L is the heat of 
vaporization for water (2.44 MJ/kg at 25 °C [24]). 


4.3. Devolatilization and char formation 

The actual burning of wood is preceded by a process called 
devolatilization. This is a series of chemical reactions, in which 
the wood decomposes into a mixture of gaseous species called 
volatiles and a solid residue called char. The combustion of 
volatiles occurs through a series of homogeneous gas-phase 
reactions, while the combustion of char occurs through 
heterogeneous oxidation reactions at the surface of the char. 

The chemical composition of the volatiles obviously differs 
from that of char. Hydrogen and oxygen in the wood tend to be 
released with the volatiles, whereas char is rich in carbon. 
This is illustrated in Table 7 which is based on the data pre¬ 
sented by Roberts [25] for European beech during different 
stages of devolatilization. It can be seen that for the volatiles, 
the values of both m and n are approximately unity; for char, 
the values of both m and n approach zero as the devolatiliza¬ 
tion proceeds. However, the value of m - n only varies in the 
range 0.03 ... 0.1 both for the volatiles and for the char. 

It has also been observed that the heating value of the 
volatiles is lower than the heating value of char [2,3], However, 
changes in heating value are not really important for deciding 
whether the oxygen consumption method and the carbon 
dioxide generation method can be applied in the analysis of 
wood firing; what really matters is whether the values of Q£ F 
and Q^'p are the same for the burning of the volatiles as they 
are for the burning of the char. The data presented by Roberts 
[3,25] seem to support the conclusion that this indeed is the 
case; see Table 8. Both for the volatiles and for the char, the 
values of Q£ F and the values of Q£' F are well within the ranges 
established in Section 4.2. 


Table 8 - Lower heating values of dry wood (European 
beech) at different stages of devolatilization [3,25]. The 
numbers after the name of the wood species indicate the 
percentage of initial weight of the sample remaining at 
each stage. 


Wood species 

Heating values of dry wood 


q" 

[kj/kg] 

Qdf 

[kj/mol] 

Q£'f 

[kj/mol] 

European beech — volatiles 

15038 

450.00 

456.84 

European beech — total 

18036 

437.01 

447.17 

European beech — 89.5% 

18476 

442.64 

452.47 

European beech — 80.0% 

19105 

437.31 

446.90 

European beech — 72.4% 

19503 

436.03 

446.63 

European beech — 56.4% 

20592 

424.58 

440.47 

European beech — 48.6% 

21789 

437.69 

446.94 

European beech — 25.8% 

26542 

422.57 

444.22 

European beech — 17.0% 

33947 

421.59 

438.96 


4.4. Summary 

Based on the discussions presented in Sections 4.1-4.3, we 
conclude that in wood firing, the values of both Q£ F and Q£' F 
are reasonably constants regardless of wood species and 
regardless of whether we are studying the combustion of 
volatiles or the combustion of char. If we do not know the 
properties of the wood, we may safely assume that 
Q£ f = 425 kj/mol and that Q£' F = 446 kj/mol. In this case, the 
uncertainty level is expected to be approximately ±10%. If we 
know the heating value and the composition of the wood, 
then it is possible to calculate more accurate estimates by 
using the equations presented in this paper. 

The conclusion regarding the effect of wood species is 
being supported by a large body of data and can be regarded as 
rather reliable. The conclusion regarding the similarities 
between volatiles combustion and char combustion is based 
on the analysis of the results of Roberts [3,25]; it would 
certainly be worthwhile to collect and analyze more data to 
check the validity of this conclusion. 

Thus, both the oxygen consumption method and the 
carbon dioxide generation method appear to offer attractive 
possibilities for the analysis of wood firing in stoves. 


Table 7 - Chemical composition of dry wood (European beech) at different stages of devolatilization [3,25]. The numbers 
after the name of the wood species indicate the percentage of initial weight of the sample remaining at each stage. 


Wood species Chemical composition of dry wood 




<* 

[] 

W>H 

[] 

[] 

1X7 Ash 

[] 

Sum 

[] 

m 

[] 

n 

n 

m — n 

n 

m + n 

[] 

European beech — 

volatiles 

0.395 

0.069 

0.536 

0.000 

1.000 

1.048 

1.018 

0.030 

2.066 

European beech — 

total 

0.484 

0.060 

0.450 

0.006 

1.000 

0.744 

0.697 

0.046 

1.441 

European beech — 

89.5% 

0.490 

0.059 

0.443 

0.008 

1.000 

0.722 

0.678 

0.044 

1.401 

European beech — 

80.0% 

0.513 

0.057 

0.426 

0.004 

1.000 

0.667 

0.623 

0.044 

1.289 

European beech — 

72.4% 

0.524 

0.056 

0.414 

0.006 

1.000 

0.641 

0.593 

0.049 

1.234 

European beech — 

56.4% 

0.561 

0.054 

0.376 

0.009 

1.000 

0.578 

0.503 

0.075 

1.080 

European beech — 

48.6% 

0.585 

0.049 

0.359 

0.007 

1.000 

0.503 

0.460 

0.042 

0.963 

European beech — 

25.8% 

0.717 

0.040 

0.222 

0.021 

1.000 

0.335 

0.232 

0.103 

0.567 

European beech — 

17.0% 

0.928 

0.015 

0.018 

0.036 

0.997 

0.097 

0.015 

0.082 

0.112 


[ ] means that a quantity is dimensionless (e.g. a mole fraction). 
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We also conclude that the variations in the chemical 
composition of wood are quite small (bearing in mind, 
however, that we are only interested in the proportions of 
carbon, hydrogen and oxygen in wood and do not care about 
the finesses of organic chemistry). In particular, we note that 
the values of the parameters m and n are expected to behave as 
illustrated in T ables 3 and 4. If we do not know the composition 
of the wood, we may safely assume that m = 0.7 and n = 0.6. 
Such an assumption will not add significantly to the uncer¬ 
tainty level of our analysis, as will be discussed in Sections 7 
and 8. Of course, if we do know the chemical composition of 
the wood, we may use accurate values for m and n. 


5. Experimental arrangement for the 
application of the methods 

In this paper, we consider the experimental arrangement 
shown in Fig. 1. Flue gas from burning fuel is collected into an 
exhaust duct, and flow rate of the flue gas and the mole 
fractions of oxygen and carbon dioxide in the flue gas are 
measured and recorded as functions of time. This information 
is used to determine the number of moles of oxygen 
consumed in the combustion reactions and/or the number of 
moles of carbon dioxide generated in the reactions. 

In the field of fire safety engineering, the experimental 
arrangement shown in Fig. 1 is called an open system. This is 
the commonly used arrangement. It would also be possible to 
build a closed system, in which the rate of air inflow into the 
system is measured rather than the rate of flue gas outflow 
from the system. However, open systems have the advantage 
that the amount and composition of the flue gas are measured 
at the same location, and air leaks therefore have no effect on 
the performance of the system. To reach a similar level of 
accuracy, a closed system would have to be essentially gas- 



Fig. 1 - Schematic view of the experimental setup for 
analysis of combustion processes using the oxygen 
consumption method and the carbon dioxide generation 
method. 


tight, which is quite difficult to achieve in practice. Thus, we 
will limit our analysis to open systems. 


6. Chemical reaction analysis 

6.1. Reactants and products 

As described in Section 4, the chemical formula of wood is 
assumed to be (CH 2m On)y, where the values of m and n are 
expected to be slightly less than unity. 

The combustion air entering into the system is assumed to 
consist of dry air (mole fraction k) and water vapor (mole 
fraction 1 - fe). These can be assumed to be known; the value 
of fe can be computed from the values of ambient pressure, 
ambient temperature and relative humidity by using, e.g., the 
procedure presented in Refs. [14,26], Normally, fe = 0.98...0.99. 
Dry air is further divided into oxygen, nitrogen and carbon 
dioxide (here nitrogen also includes argon and other inert 
gases that are not of interest for our purposes). The mole 
fractions of oxygen, nitrogen and carbon dioxide in dry air are 
denoted by/, g and h, respectively. In numerical calculations 
we will use the values/= 0.2095 and h = 380 lCT 6 , but this 
does not limit the generality of our approach, and any other 
values could indeed be used if necessary. 

We also assume that the combustion is complete, that 
is, the wood is oxidized to carbon dioxide and water vapor 
and that no soot, carbon monoxide, unburned hydrocar¬ 
bons, or nitrogen oxides exist in the flue gas. Thus, the flue 
gas will consist of carbon dioxide, water vapor, oxygen and 
nitrogen. 

6.2. Reaction equation and solution strategy 

Let us ignore the time-dependent nature of the burning 
process for a while. Thus, the chemical reaction equation for 
the burning of dry wood in air can be written as 

(CH 2m O„) y + Z{fe[f0 2 + gN 2 + hC0 2 ] + (1 - fe)H 2 0} 

->pC0 2 + sH 2 0 +10 2 + uN 2 

There are a total of four unknowns in Equation (13): m, n, y 
and Z. The first three of them specify the amount and chem¬ 
ical composition of wood participating in the reaction, and Z 
specifies the amount of combustion air. At this point, the 
stoichiometric coefficients of the products (p, s, t and u) are 
also unknown, but they can be expressed in terms of input of 
the reactants, as shown in Table 9. 

So we need to solve the values of four unknowns, but there 
are only three measurements in the experimental arrange¬ 
ment shown in Fig. 1. Thus, there will be only three equations 
available for solving the values of the four unknowns. Some 
uncertainty will therefore remain in the final outcome of the 
analysis. This is normally accepted, but it is worth noting that 
the uncertainty could actually be eliminated by adding one 
more measurement to the system. For example, the 
measurement of the water vapor mole fraction in the flue gas 
has been advocated in Refs. [12,13,26], Another option is to 
measure the mole fraction of argon in flue gas and to compare 
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Table 9 - 

Elementbalance for left and right-hand sides of Equation (13) and solution of stoichiometric coefficients p, s, t and u. 

Element 

Reactants 

Products 

Solution 

C 

y + Zkh 

P 

p = y + Zfeh 

h 2 

y m + Z(1 — fe) 

s 

1 

s = ym + Z(1 — fe) 

o 2 

^yn + Zkf + Zkh + ^Z( 1 - fe) 

P + 2 S + t 

t = Zfe/-y| 1 + i(m-n)| 

n 2 

Zkg 

u 

u = Zkg 


that with the mole fraction of argon in combustion air (which 
can normally be assumed to be known) [6]. 

We will now concentrate our attention to the problem with 
four unknowns and three equations. Intuitively, one would 
assume that y and Z must be eliminated from the solution 
(since there is no way to estimate their values) but that m and 
n can be retained (since we have seen that we do know quite 
a few things about their values). This is the strategy that will 
be adopted in our study. 

6.3. The amount and composition of flue gas 

The flow rate measured in the test arrangement shown in 
Fig. 1 is the flow rate of the wet flue gas (which includes water 
vapor). On the other hand, the mole fractions of 0 2 and C0 2 
are usually measured on a dry basis, that is, water vapor is 
removed from the flue gas before it flows to the gas analyzers. 
Thus, it will be necessary to distinguish between wet flue gas 
(which includes water vapor) and dry flue gas (from which 
water vapor has been removed). This is achieved by using 
subscripts WFG and DFG. 

The first step is to derive expressions for the amount of flue 
gas (both dry and wet will be needed) and for the mole frac¬ 
tions of oxygen and carbon dioxide in dry flue gas. By using 
Table 9, the number of moles of dry flue gas is 


N D fg = p +1 + u 

= y + Zfeh + Zkf - y 


l--(m-n) 


+ Zkg 


= Zk- —y(m - n) 


and the number of moles of wet flue gas is 


(14) 


Nwfg — Nqfg + S 


= Zfe — -y(m — n) + y m + Z(1 — fe). 


(15) 


= Z + -y (m + n) 

The mole fraction of 0 2 in dry flue gas is 


zkf- y 


Xo, = 


1 + 2 ( m — ”) 


Nn 


Zfe - -y(m - n) 


(16) 


6.4. Consumption of oxygen in the reaction and 
estimation of the heat release rate 


We can now write the expression for the oxygen consumption 
in the reactions: 


ANo 2 — No 2 out — No 2 ,in 

= Xq 2 N d fG — Zkf 


= Xo, 


Zk - -y (m - n) 


-Zkf 


= ZkAx 0 , - ^y(m - n)x 0 . 


where 


( 18 ) 


Ax 02 =x 0 ,-f (19) 

is the difference between the oxygen mole fraction in dry flue 
gas and oxygen mole fraction in dry combustion air. 

The next step is to eliminate y from Equation (18). Equation 
(16) gives 


y = 


-ZkAx 0 


l + -(m-n)(l-x 02 ) 


( 20 ) 


and substitution into Equation (18) gives 
1 -ZfeAxo 


l + ^(m-n)(l —Xq 2 


(m - n)x 02 


= ZfeAxo 


l + -{m-n) 
l + ^(m-n)(l —Xo 2 


( 21 ) 


We continue by eliminating Z. Equations (15) and (20) are 
used to obtain 


Nwfg = Z + 2 ( m + n )y 


= Z + -(m + n)-j 


-ZkAx 0 


= Z- 


l + -(m-n)(l-x 02 ) 

1 1 

l + j(m-n)(l- Xo 2 ) --fe(m + n)Ax 02 

l+-(m-n)(l-x 02 ) 


( 22 ) 


from which we can solve 


and the mole fraction of C0 2 in dry flue gas is 


Xco 2 


p y + Zfeh 

[dfg Zfe - ^y(m - n) 


(17) 


Z = No 


1 

1 + 2 ( m 


^(I-Xq,) 


1 1 

1 + - (m - n) (1 - x 02 ) - ^k(m + n)Ax 0 , 


(23) 


and substitution into Equation (21) gives 
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1 

l + -(m-n) 

ANo 2 = N wf x;feAx 0; -jj- 

l + -(m-n)(l-Xo 2 ) --k(m + n)A xo 2 

(24) 

Our main interest, however, is in non-steady systems where 
the flue gas flow rate and composition change with time. 
Thus, we express the instantaneous rate of oxygen 
consumption with the equation 


1 

l + -(m-n) 

ANo 2 = N wrG feAXo ; - jj- 

l + -(m-n)( 1 -Xo 2 ) --fe(m + n)Ax 02 

(25) 

where Nwfg is the instantaneous molar flow rate of wet flue 
gas in the exhaust duct. The measured value of the oxygen 
mole fraction is also to be understood as an instantaneous 
value. 

Thus, we are able obtain the instantaneous value of the 
heat release rate as 


Q=-ANo 2 Q" 
= —NwFcfeQ" ' 


1 + 2 (m-n) 


Ax, 


o 2 


(26) 


1 1 

1 + ;r(m - n)(l - x 02 ) --fe(m + n)Axo 2 


Here a negative sign was introduced in order to make Q 
positive (ANo 2 is negative and Q" is positive). 

Bearing in mind that we are discussing the combustion of 
dry wood, we may note that Q" = Q£ F = 425 kj/mol and 
rewrite Equation (26) as 


Q — —NwrcfeOcF 


l + -(m- n ) 


Ax 0 


1 1 

l + -(m-n)(l-Xo 2 ) - -k(m + n)A x 02 


(27) 


Equation (27) is the basic equation for the oxygen 
consumption method. Essentially, what we have accom¬ 
plished so far is that we first formulated two equations using 
the data from two measurements (x 02 and N WF g)- The equa¬ 
tions were then used to eliminate two of the unknowns in our 
original problem (y and Z). Two more unknowns are still 
remaining and must somehow be dealt with: m and n. For this 
purpose, we have one more measurement (x C o 2 ), which can be 
used to say something about the values of m and n but not to 
solve the values of the both of them. However, before we get 
into that we note that Equation (27) is actually quite powerful 
even without any further manipulations: if we are able to 
estimate the values of m — n and m + n, then we can use 
Equation (27) to solve for the heat release rate, and the 
measurement of the mole fraction of C0 2 in flue gas is not 
even needed. 

In the case of wood firing, the values of m - n and m + n can 
be estimated based on the information presented in Section 4. 
The value of m - n was found to be approximately 0.1 both 
during the burning of the volatiles and during the burning of 
the char, but the value of m + n was found to range from 
approximately 0 (during char burning) to approximately 2 
(during volatiles burning). Some uncertainty will therefore 


remain in the results. The possibility of decreasing the 
uncertainty by the measurement of the carbon dioxide mole 
fraction will be discussed in Section 6.6. 

The best results for wood firing analysis can probably be 
obtained by setting m + n = 1.3. This is the average value for all 
wood species analyzed in Section 4. Using the average value 
for wood means that the consumption of oxygen (and thus the 
heat release rate) will be overestimated during the burning of 
the volatiles (since the value of the denominator in Equation 
(27) will then be smaller than the actual value) and under¬ 
estimated during the burning of the char (since the value of 
the denominator in Equation (27) will then be larger than the 
actual value). The total amount of energy released during 
combustion is obtained by time integration of heat release 
rate over the whole duration of burning and can be expected 
to be close to the actual value. 

If the average value of m + n is used in the estimation of the 
heat release rate rather than the actual value of m + n, an error 
will be introduced. The magnitude of the error can be estimated 
as follows. Consider the char burning phase and assume 
that xo 2 =0.1 (as would be expected for good stoves). Thus, 
Axo 2 = - 0.1 and 1 - x 02 =0.9; these values are substituted into 
Equation (27). The first term in the denominator in the right- 
hand side of Equation (27) is unity. The second term can be 
expected to be approximately 0.045, since m-n = 0.1. The 
actual value of the third term is zero, (since m + n = 0 during 
char burning) but the value used in the computations will be 
approximately 0.5 ■ 0.99 • 1.3 ■ (-0.1) = -0.064. The actual value 
of the denominator is therefore 1 + 0.045 = 1.045, whereas the 
value used in the computations is 1 + 0.045 - (-0.064) = 1.11. 
The relative error is in this case approximately -6%. 

In the next section, we will formulate an analogous 
expression for computing the heat release rate based on the 
number of moles of carbon dioxide generated during the 
burning. After that, we will get back to the refinement of 
Equation (27) in Section 6.6, and, finally, analyze the effect of 
fuel moisture in Section 7. 


6.5. Generation of carbon dioxide in the reaction and 
estimation of the heat release rate 


Similarly to what was carried out in the previous section, we 
write the expression for the carbon dioxide generation in the 
reactions: 


ANco 2 — Nco 2 ,out — N C o 2l in 
= Xco 2 N DF g — Zfeh 


— Xco 2 


Zfe - -y (m - n) 


Zfeh 


= ZfeAx C o 2 - ^y(m - n)x C o 2 


(28) 


where 


Ax C o 2 = x C o 2 - h (29) 

is the difference between the carbon dioxide mole fraction in 
dry flue gas and carbon dioxide mole fraction in dry 
combustion air. 

The next step is to eliminate y from Equation (28). Equation 
(17) gives 
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y = - 


ZfeAXc 


1 + 2 (m - n)x C o 2 

and substitution into Equation (28) gives 


(30) 


ANco 2 — ZfeAXco, — — ■ 


ZfeAXr 


1 + j ( m - n ) x co 2 


- (m - n)x C o 2 


(31) 


= ZfeAXr 


1 + 2 (m - n)x C o 2 


We continue by eliminating Z. Equations (15) and (30) are 
used to obtain 


Nwfg = Z + 2 ( m + n )y 

_ . 1 , , , ZfeAxco. 

= Z + ^ (m + n)-j-- 


= Z- 


1 + j ( m - n ) x co 2 

1 1 

1 + - (m - n)x C02 + 7 rfe(m + n)Ax C02 


(32) 


1 + j ( m - n ) x co 2 


from which we can solve 


1 

l + -(m-n)x C02 

z = Nwfg -j-^-j- (33) 

1 + - (m - n)x C02 +-k(m + n)Ax COl 

and substitution into Equation (31) gives 

1 

AN C o 2 = NwFGfeAx C o 2 -j-j- (34) 

1 + -(m - n)x C02 + -fe(m + n)Ax C o 2 

Again, our main interest is in non-steady systems where the 
flue gas flow rate and composition change with time. Thus, we 
express the instantaneous rate of carbon dioxide generation 
with the equation 


1 

AN(;o 2 = N wrG feAx C o 2 -j-j- (35) 

1 + 2 ( m - tl)x C o 2 + 2 te ( m + n )Ax C o 2 

where Nwfg is the instantaneous molar flow rate of wet flue 
gas in the exhaust duct. The measured value of the carbon 
dioxide mole fraction is also to be understood as an instan¬ 
taneous value. 

The instantaneous rate of heat release is therefore ob¬ 
tained as 


Q = AN c02 -Q'" 

= NwFGkQ"'-j-- (36) 

1 + 2 ( m - n)x C o 2 + 2 fe ( m + Ax co 2 

and, since we are discussing the combustion of dry wood, we 
note that Q'" = Q£' F = 446 kj/mol (as was shown earlier). Thus, 
Equation (36) may be rewritten as 

A y 

Q = NwFGfeQ^'F-J-^- (37) 

1 + 2 ( m - n ) x co 2 + 2 fe ( m + n )Ax C o 2 

Equation (37) is the basic equation for carbon dioxide 
generation method. Analogously to what was done during the 


derivation of Equation (27), we first formulated two equations 
using the data from two measurements (this time they were 
x C o 2 and Nwfg)- The equations were then used to eliminate 
two of the unknowns in our original problem (y and Z). What 
we have achieved here is quite similar to what was achieved 
in the previous section. Again, we have two more unknowns 
to be dealt with (m and n) and one more measurement (this 
time Xo 2 ), which can be used to say something about the 
values of m and n but not to solve the values of the both of 
them. 

Also we note that Equation (37) is quite powerful even 
without any further manipulations: if we are able to estimate 
the values of m~ n and m + n, then we can use Equation (37) to 
solve for the heat release rate and the measurement of the 
mole fraction of 0 2 in flue gas is not even needed. 

As was noted in the previous section, the value of m - n is 
approximately 0.1 both during the burning of the volatiles and 
during the burning of the char, but the value ofm + n varies 
between approximately 0 (during char burning) and approxi¬ 
mately 2 (during volatiles burning). Some uncertainty will 
therefore remain in the results. The possibility of decreasing 
the uncertainty by the measurement of the oxygen mole 
fraction will be discussed in Section 6.6. 

The best results for wood firing analysis can probably be 
obtained by setting m-n = 1.3. This is the average value for all 
wood species analyzed in this paper. Using the average value 
for wood means that the generation of carbon dioxide (and 
thus the heat release rate) will be overestimated during the 
burning of the volatiles (since the value of the denominator in 
Equation (37) will then be smaller than the actual value) and 
underestimated during the burning of the char (since the 
value of the denominator in Equation (37) will then be larger 
than the actual value). The total amount of energy released 
during combustion is obtained by time integration of heat 
release rate over the whole duration of burning and can be 
expected to be close to the actual value. 

If the average value of m + n is used in the estimation of the 
heat release rate rather than the actual value of m + n, an error 
will be introduced. The magnitude of the error can be esti¬ 
mated as follows. Consider the char burning phase and 
assume that x C o 2 =0.1 (as would be expected for good stoves). 
Thus, Ax C o 2 =0.1. The first term in the denominator in the 
right-hand side of Equation (25) is unity. The second term can 
be expected to be approximately 0.005, since m - n~ 0.1. The 
actual value of the third term is zero (since m + n = 0 during 
char burning), but the value used in the computations will be 
approximately 0.5 ■ 0.99 • 1.3 • 0.1 = 0.064. The actual value of 
the denominator is therefore 1 + 0.005 = 1.005, whereas the 
value used in the computations is 1 + 0.005 + 0.064 = 1.07. The 
relative error is in this case approximately -6%. 

In the next section, we will analyze what can be achieved 
when both the oxygen mole fraction and the carbon dioxide 
mole fraction are used in the analysis. 

Finally, it may be noted that Equation (35) can be used to 
compute the total amount of carbon dioxide generated during 
the test. If the quantity and composition of the fuel are 
known, this information can be used as an additional check 
on the reliability of the results. The equations for computing 
the total amount of carbon dioxide generated during the 
test are 
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AN C o 2 — J AN C o 2 dt (38) 

for the case where measurement signals are continuously 
recorded (as might be the case when an analogous system is 
used for data logging) and 

AN C o 2 = N AN c q 2 At (39) 

for the case where the measured values are discretely recor¬ 
ded using a time interval At (this would normally be the case 
when a digital system is used for data logging). 


6.6. Analysis of fuel composition from measurements of 

flue gas composition 


In the previous two sections, we derived the basic equations 
for the oxygen consumption method and for the carbon 
dioxide generation method. In each case, it was assumed that 
the experimental data that are available consist of the total 
amount of flue gas and the mole fraction of one component of 
the flue gas: in the case of the oxygen consumption method 
this component was oxygen, and in the case of the carbon 
dioxide generation method this component was carbon 
dioxide. Normally, however, the experimental data include 
both the mole fraction of oxygen and the mole fraction of 
carbon dioxide. In this section, we shall see what improve¬ 
ments this fact can be expected to bring. 

The problem that we have in our hands is that we do not 
know the values of the parameters m and n. Of course, we 
cannot expect that the addition of one more measurement 
will make it possible to solve the values of the both of them. 
However, we have already seen the value of m - n can be 
assumed to be constant, and that it would therefore be suffi¬ 
cient to be able to deduce the value of m + n from the exper¬ 
imental data. In the following we shall see whether this is 
possible or not. 

The simplest way to proceed is to note that the difference 
between the mole fraction of oxygen in dry flue gas and the 
mole fraction of oxygen in dry air is 


Axq 2 


*o 2 -/ = 


y 


l + i(m-n)(l-/) 

1 I 1 

Zk - -y(m - n) 


(40) 


and the difference between the mole fraction of carbon 
dioxide in dry flue gas and the mole fraction of carbon dioxide 
in dry air is 


Ax C o 2 = x C q 2 - h = - 


1 + -(m - n)h 


Zk - -y(m - n) 


(41) 


Now we can define a new parameter £ with the equation 


the combustion air, not on the amounts of fuel and combus¬ 
tion air participating in the reaction. 

Rearrangement of Equation (42) gives 


m 


2 ( 1 - 

«i -/) - 


2 1 + 


Ax C0 A 
Axq 2 ) 


Ax C q 2 

AXq 2 


(1 ~f)~h 


(43) 


which can be inserted into Equations (27) and (37) to eliminate 
m - n. It is obvious that both the oxygen consumption method 
and the carbon dioxide generation method will now converge 
into a single method. This is due to the fact that although the 
routes are different, the end results are based on exactly the 
same information, that is, on the measured values of the flue 
gas flow rate and the mole fractions of oxygen and carbon 
dioxide in the flue gas. Thus, if we start from Equation (27), we 
will obtain the heat release rate as a function of Nwfg, k, Qjj F , 
x 0! , x co 2 and m + n; if we start from Equation (37), we will 
obtain the heat release rate as a function of Nwfg, k, Qjj' F , x 02 , 
x C o 2 and m + n. These two expressions can be shown to be 
identical when we observe that there exists a relationship 
between Q£ F and Qjjj. (see Section 4). We do not go through the 
arithmetics here, but in the next section we will present an 
alternative approach for utilizing the measured values of both 
Xo 2 and x C o 2 to obtain an expression for the heat release rate. 
The starting point will be the basic equation of the 
oxygen consumption method, and the result will be an 
equation quite similar to what has been used in the field of fire 
safety science. 

We have now come as far as we can regarding the possi¬ 
bility of using the experimental data that we possess to 
determine the composition of the fuel. In other words: the 
addition of another measurement of flue gas composition has 
not really helped us at all. What we get from the additional 
measurement is the value of m - n, which we already know 
(based on what was presented in Section 4); what we do not 
get is the value of m + n, which we hoped to be able to 
determine. 

How and why did things go this way? 

By looking at Equations (16) and (17), we note that m and 
n do not appear independently from each other; they only 
appear as the expression m - n. Thus, no matter what 
manipulations are carried out, we will only be able to 
solve the value of m — n. What was shown in Equations 
(40)—(43) is just the simplest and most elegant way to 
accomplish this. 

It is also easy to see what will be needed to achieve more 
useful results: we will need an equation where either m or n 
appears independently or where m and n appear as some 
other combination than m - n. For example, consider the 
addition of the measurement of water vapor mole fraction, as 
was discussed in Section 6.2 (naturally, this is the mole frac¬ 
tion of H 2 0 in uiet flue gas). By looking at Table 9 and Equation 
(15), we obtain the equation 


Ax C 02 l + ^(m-n)h 
AX ° J l + ^(m-n)(l-_f) 


(42) 


x h 2 o = I 


ym + Z(1 - fe) 
1 

Z + -y(m + n) 


(44) 


Thus, we have been able to eliminate y and Z. The value of £ 
will therefore only depend on the compositions of the fuel and 


which can be used to solve for m + n as a function of m. 
Adding together the expression for m + n and the expression 
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for m - n from Equation (43) will result in an equation from 
which m can be solved. The arithmetics that are needed will be 
considerable and are not shown here. 

In practice, the following options are available: (1) measure 
the mole fractions of oxygen and carbon dioxide in dry flue gas 
and the mole fraction of water vapor in wet flue gas, and solve 
both m and n from the measurement data; (2) measure the 
mole fraction of oxygen in dry flue gas and the mole fraction of 
water vapor in wet flue gas, set m - n = 0.1 based on what was 
presented in Section 4 and then solve m + n from the 
measurement data; or (3) measure the mole fraction of carbon 
dioxide in dry flue gas and the mole fraction of water vapor in 
wet flue gas, set m - n = 0.1 based on what was presented in 
Section 4 and then solve m + n from the measurement data. 
Other combinations of measurements could also be 
employed. 

Here it is worth noting that the oxygen consumption 
method and the carbon dioxide method were originally 
developed for use in fire safety research. However, there are 
two important differences between experiments related to 
wood firing in stoves and experiments related to accidental 
fires. Firstly, while we know quite a lot about the composition 
of firewood, we usually do not know that much about the 
composition of the materials which burn in accidental fires; 
thus, in fire experiments it is a good practice to measure the 
mole fractions of both oxygen and carbon dioxide in flue gas. 
This will give us information of m - n, which would other¬ 
wise be quite difficult to estimate. Secondly, the values of 
Axo 2 and Ax C o 2 are quite different in the flue gas from wood 
firing in stoves and in the exhaust gas from fire experiments. 
In the case of wood firing in a stove, we generally try to 
minimize the amount of excess air in order to maximize the 
energy efficiency of the stove; thus, the absolute values of 
Axo 2 and Axco 2 are expected to be of the order of 0.1. In the 
case of fire experiments, the combustion products are 
mixed with a very large amount of air before they enter the 
exhaust duct; thus, the absolute values of Ax 02 and Ax C o 2 are 
often below 0.01 throughout the whole experiment and 
seldom larger than 0.03. Consequently, in fire experiments 
the uncertainty caused by the fact that we do not know 
and cannot even determine the value of m + n is much 
less significant than in experiments related to wood firing 
in stoves. 

As a final note for this section, we point out that the value 
of £ can be used to check the operation of the gas analyzers if 
the composition of the fuel is known. This possibility is also 
noted in [14]. The idea is that if the value of £ computed from 
the experimental data differs from the theoretical value 
computed from the fuel composition, then something must be 
wrong. The accuracy of modern gas analyzers tends to be 
excellent, and this information can indeed be used for diag¬ 
nostic purposes. 

The values of £ for the organic compounds considered in 
this paper are given in Table 10. It can be seen that for these 
compounds, the value of £ only varies between approxi¬ 
mately 0.56 and 1. The value of £ for propane is 0.655; for 
some reason, an erroneous value 0.6 is given in [14]. For dry 
wood, m - n ~ 0.1 which yields £ = 0.96. We will use this 
information later in Section 9 for checking the validity of our 
gas analyses. 


Table 10 - The values of parameter £ for the organic 
compounds considered in this paper. The computations 
were carried out with / = 0.2095 and h = 380 10“ 6 . 

Fuel 

Chemical formula 

m 

n 

£ 

Carbon 

C 

0 

0 

1 

Methane 

ch 4 

2 

0 

0.5587 

Propane 

c 3 h 8 

4 

3 

0 

0.6550 

Benzene 

C 6 H 6 

1 

0 

0.8351 

Methanol 

ch 3 oh 

2 

2 

1 

0.7169 

Acetic acid 

ch 3 cooh 

1 

1 

1 

D-glucose 

c 6 h 12 o 6 

1 

1 

1 


6.7. Oxygen consumption method and the formalism 
employed in fire safety studies 

In the field of fire safety engineering, the application of the 
oxygen consumption method is normally based on the equa¬ 
tions presented by Parker [11,12]. The outward appearance of 
those equations is quite different from that of Equation (27), 
which was derived in the Section 6.4. Here we start from 
Equation (27) and derive an expression which allows one to 
see the connection between our approach and Parker’s 
equations. The algebraic manipulations are actually quite 
remarkable and are presented here in a shortened form. 

We start by writing 



1 




AXq 2 


1 + 2 ( m ~ n ) (1 - x o 2 ) - 2 fe ( m + n )Axo 2 


and define an auxiliary parameter 

Q = 1 — h — £/ (46) 

where £ is the ratio of Ax C o 2 to -Ax 02 as was defined in Equa¬ 
tion (42). After some algebra, we obtain 


(1-/-K) 


£2 = 


l + -(m-n) 


l + o(m - n)(l -/) 


and 


(47) 


£2Ax 02 = -[/(l - x C o 2 ) - Xo 2 (l - h)] 

The next step is to write 


R = R- 

where 

Q T 

1 


S = Q 

1 + 2 ( m “ n ) 

Ax 02 

= - 

1 

1 + 2 ( m - n) 

' [f(l _ x co 2 ) 

and 



T = r + T" 



(48) 


(49) 


(50) 


( 51 ) 
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where 


T = Q- l + -(m-n)(l-x 02 ) (52) 

and 

T" = Q- ~^(m + n)feAxo 2 (53) 

To solve for T, we note that it can be shown that 


l + -(m-n)(l-Xo 2 ) = 


l + -(m-n)(l-/) 


(l ~ Xq 2 — X C q 2 ) 


1-f-h 


1 + 2 ( m — n ) 


(1 — Xo 2 — x C o 2 ) 


and substitution into Equation (52) gives 


T = 


l + -(m-n) 


(l — x ( ) 2 — Xco 2 ) 


(54) 


(55) 


To solve for T", we note that Equation (53) can be also 
written as 


T" = QAx a , 


-(m + n)k 


(56) 


= [/(l-x C o 2 ) -Xo 2 (l-h)] --(m + n)fe 
Now we are ready to reassemble Equation (49). We obtain 


- r w 

-(m + n) 

1 + fe • ^ '/0 

l+ 2 ( m “ n ) 

(59) 

Secondly, the fuel composition is 
parameters a and /3. For a fuel having 
C a H b O c , these parameters are defined by 

expressed by using 
a chemical formula 

4a + 2b 

4a + b - 2c 

(60) 

and 


a = 1 + kf(/3 — 1) 

(61) 

In our case, a = y, b = 2my and c = ny. 

Thus, 

f> = -- 

l + -(ni-n) 

(62) 

and 


1. 

-(m + n) 
a = l + kf 2 1 

l + -(m-n) 

(63) 

Substitution into Equation (59) gives 


P _ -/0 

1 + (a - 1)0 

(64) 


and Equation (27) becomes 


S _ S 
T - T + T" 


- 

1 

1 + -(m-n) 

■ [f (1 — x co 2 ) — x o 2 (1 — h)] 

1 

1 + - (m - n) 

(l — Xo 2 — Xc 0 2 ) + [f (l — Xco 2 ) — Xo 2 (1 — h)] • 

1 

- (m + n)fe 


/(l-Xco 2 ) -Xp 2 (l-h) 

1 — Xo 2 — Xco 2 

1|h ^ (m + n) /(1 —x C02 )-x 02 (l-h) 

l+i( m _„) 1-Xo 2 -Xco 2 


(57) 


In the field of fire safety engineering, the notation is 
somewhat different from the notation adopted in this study. 
The following two steps are therefore needed. 

Firstly, the composition of the flue gas is usually expressed 
in terms of the oxygen depletion factor, which is defined by 


/(1 — Xco 2 ) — Xo 2 (l — h) 

f( 1 — Xq 2 — x C o 2 ) 


(58) 


This can be obtained from Equation (28) of Parker [11] by 
setting the mole fraction of carbon monoxide equal to zero 
and by assuming that the gas flowing to the gas analyzers is 
dry. Substitution of Equation (58) into Equation (57) gives 


Q = -NwFGkQ£ F R 


-J0 


- - NwFGfeQ ^TTF"-1)0 

0 


— Nwfg H/Qdi 


1 + (a - 1)0 


(65) 


We have now arrived at something that resembles the 
equations employed in the field of fire safety engineering; see, 
e.g., [14,13]. The main difference is that in fire safety engi¬ 
neering it has become a common practice to express flow 
rates as volume flow rates rather than molar flow rates. This is 
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in accordance with the fact that the heating value is also 
expressed in terms of the volume of oxygen consumed in the 
combustion reactions, not in terms of the number of moles of 
oxygen consumed in the combustion reactions. This point 
was already noted earlier. The use of volume flow rates can be 
criticized, since it requires that the reference temperature and 
pressure must be specified, which provides a potential source 
of confusion. This point was noted by Janssens [26], who chose 
to use mass flow rates in computations. We have chosen to 
use molar flow rates, which in our opinion give the simplest 
way to express the connection between flue gas analysis and 
oxygen consumption. 

Another point to be noted is that in Equation (65), every¬ 
thing else can be measured or is known with good accuracy 
except the value of a, which depends on the composition of 
the fuel. In the following, we will analyze the magnitude of 
this uncertainty. 

Let us first separate the effect of air humidity from the 
effect of fuel composition. Define 

1, , 

-(m + n) 

a 1 = 1 +f{P — 1) = 1 +_f —=-j- (66) 

l + -(m-n) 

Substitution of Equation (66) into Equation (59) gives 


R_ -/* 

1 + fe(a' — 1 )</> 

and Equation (27) becomes 


(67) 


Q.DF — “NwFGfeQ() F R 

= “N WFG feQ£ ri + fe( ^_ 1)0 (68) 

= NwPGfefQS Fl + fe( J _ 1)0 

Let us now take a look at the term k(o:' - 1 )</> in the denomi¬ 
nator of Equation (68). The value of k will be approximately 
unity and the value of tp is expected to be approximately 0.5 for 
good stoves but not higher than approximately 0.2 during fire 
experiments. The values of a' for the organic compounds 
considered in Section 2.1 are given in Table 11. It can be seen 
that all values are within the range from 1.0 to 1.21. For the 
wood species considered here, the average value of m + n was 
approximately 1.3 and the average value of m - n was approx¬ 
imately 0.1; thus, a/ = 1.13. However, it was also observed that 


Table 11 - The values of parameter a 
compounds considered in this paper, 
were carried out with / = 0.2095. 

for the organic 

The computations 

Fuel 

Chemical formula 

m 

n 

a' 

Carbon 

C 

0 

0 

1.000 

Methane 

ch 4 

2 

0 

1.105 

Propane 

C 3 H s 

4 

3 

0 

1.084 

Benzene 

c 6 h 6 

1 

0 

1.042 

Methanol 

ch 3 oh 

2 

2 

1 

1.210 

Acetic acid 

ch 3 cooh 

1 

1 

1.210 

D-glucose 

c 6 H 12 o 6 

1 

1 

1.210 


the value of m + n actually varies between 0 (during char 
combustion) and 2 (during volatiles combustion), whereas 
m - n ~ 0.1 during the whole process. Thus, in wood combus¬ 
tion the value of a / varies within the range from 1.0 (during char 
combustion) to 1.2 (during volatiles combustion). 

Thus, the value of k (a 1 - 1 )<p is expected to be within the 
range 0...0.1 for wood firing in stoves and within the range 
0...0.04 for fire experiments. By setting a' = 1.1, the uncer¬ 
tainty in the denominator of Equation (68) will be of the order 
of ±5% for wood firing in stoves and of the order of ±2% for fire 
experiments. This is the uncertainty caused by the fact that 
our measurements do not contain a sufficient amount of 
information to fully specify the composition of the fuel. A 
magnitude of the order of ±5% for the uncertainty was already 
noted in Sections 6.4 and 6.5. The fact that higher uncertainty 
must be accepted when results from experiments related to 
wood firing in stoves are analyzed than when results from 
experiments related to accidental fires are analyzed was also 
noted in Sections 6.4 and 6.5. 

Finally, it may be noted that the dependence of a on air 
humidity was apparently ignored when the oxygen 
consumption method was adopted into international fire test 
standards. In fire tests, the composition of the fuel is usually 
not known, and an average value must be chosen for a. The 
average value chosen in [14,13] is 1.105. This value is in the 
middle of the range of tabulated values of a presented by 
Parker [12,11] and was also recommended in [12,11], However, 
the tabulated values of a were for fuels burning in dry air, and 
somewhat smaller values should be used for fuels burning in 
air containing moisture. In any case, the effect of ignoring the 
dependence of a on the humidity of combustion air is negli¬ 
gible; an error of a few percent in the value of ( a - 1 )<p will not 
be important considering the fact that the value of (a - l)<j> is 
always small when compared with unity. 


7. Effect of moisture in wood 

So far, we have assumed the firewood to be completely dry. In 
this section, we shall see how moisture in the wood affects the 
determination of the heat release rate. This discussion has 
been postponed until now since we will have to deal with the 
effect of moisture both on the heating value of the wood and 
on the amount and composition of flue gas from the burning 
of the wood. The effect on heating value will be discussed in 
Section 7.1 and the effect on the amount and composition of 
flue gas will be discussed in Section 7.2. 

Note that the moisture content of the fuel is here defined as 

w , = rn h^o = _ n^o _ (69) 

Hz ° m F m H2 o + m DF 

where the subscript H 2 0 refers to moisture, F to fuel (which 
includes dry fuel and moisture) and DF to dry fuel. Throughout 
this paper, the superscript * is used to indicate quantities that 
are expressed relative to the mass of fuel. 

7.1. Effect on heating values 

This section contains an analysis of the effect of fuel moisture 
on Q" and Q'"; here Q" is the heatingvalue of the fuel per 1 mol of 
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oxygen consumed in the reaction and Q'" is the heating value of 
the fuel per 1 mol of carbon dioxide generated in the reaction. 
The heating value of the fuel is 


— (1 u Jh 2 o) ( 1df tu H20 L 
w 


= ( 1 - 


"h 2 o 


)<Jd 


h 2 o 


) <1df 


(70) 


where q DF is the heating value of the dry fuel and L is the heat 
of vaporization of water. The value of L is 2.44 MJ/kmol at 25 °C 
(which is normally the reference temperature in the deter¬ 
mination of the heating value of fuels [24]). As the values of q DF 
are typically around 19 MJ/kg (see Section 4), the value of L/q DF 
is expected to be approximately 0.13, and it can be seen that 
the latter term in the brackets in Equation (70) is quite small 
when compared with the first term except at very high values 
of moisture content. 

Analogously to the derivation of Equations (8)-(ll), we 
obtain the following results. The number of moles of oxygen 
consumed during the complete combustion of 1 kg of fuel is 



Fig. 2 - The effect of fuel moisture on Qf and Qf' for the 
hypothetical wood described in Section 7.1. Here the lower 
heating values are considered. [ ] means that a quantity is 
dimensionless (e.g. a mole fraction). 


AN* 2 = (1-u;* 20 ){aN"} (71) 

and the amount of heat released during the combustion 
divided by the number of moles of oxygen consumed in the 
reaction is 


Of = 


anj 2 

<1df 


AN, 


o. 


1- 


u H 2 0 


1 ” ^HjO <1 dF. 


(72) 


= Q£f 


w h 2 o 


1 W ihO ^DF 


The number of moles of carbon dioxide generated during 
the complete combustion of 1 kg wood is 


(CH 2m O„) y + T H 2 0 + Z{fe[f0 2 + gN 2 + hC0 2 ] + (1 - fe)H 2 0} 
—>pC0 2 + sH 2 0 +10 2 + uN 2 

which can also be written as 


(75) 


(cH 2meff 0„ eff ) y + Z{fe[f0 2 + gN 2 + hC0 2 ] + (1 - fe)H 2 0} 
->pC0 2 + sH 2 0 + t0 2 + uN 2 

where 


(76) 


m eff = m + - 


(77) 


and 


ANco 2 — (1 ^hiiO^Nco,, (73) 

and the amount of heat released during the combustion 
divided by the number of moles of carbon dioxide generated in 
the reaction is 


Qf' = 


?F 


ANco 2 

qDF 


AN; 


co 2 


= Q£'f 


1 -- 


1 ~1 

w 


w h 2 o 


) ^DF. 


h 2 o 


) ^DF. 


(74) 


Fig. 2 illustrates how the values of Q f and Qf change with 
the moisture content of wood. The heating values are 
computed for a hypothetical wood with Q£ F = 425 kj/mol, 
Q£' f = 446 kj/mol, m = 0.7 and n = 0.6. These values are 
representative of the wood species described in Section 4. It 
can be seen that the effect of moisture content is very weak 
until the moisture content exceeds approximately 30%. 


7.2. Effect on the amount and composition of flue gas 

To describe the combustion of wood containing moisture, 
water must be added to the reactant side of the chemical 
reaction expressed in Equation (13). Thus, we may write 


n e ff = n + - (78) 

y 

At this point, it may be noted that the evaporation of 
moisture presumably takes place in the early stages of 
combustion, whereas Equation (75) only describes the situa¬ 
tion in some average sense. A similar problem was already 
encountered when the burning of volatiles was compared 
with the burning of char (see Section 4). In the following, we 
will nevertheless assume that Equation (75) can be applied at 
each and every instant of time. 

To apply the basic equations of the oxygen consumption 
method and the carbon dioxide generation method, we only 
need m - n and m + n. From Equations (77) and (78), we can 
solve 

m eff - n eff = m - n (79) 

and 


ITleff + tl e ff 


m + n + 2- 

y 


(80) 


Thus, we come back to the problem that was already 
discovered in Section 6.6. The moisture in the fuel has no 
effect onm-n (the value of which can be determined from 
flue gas composition) but increases the uncertainty in the 
value ofm + n (which we cannot determine from our data). 
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To be able to continue, we need to solve the value of y/y. 
Proceed in the following manner. Let m F be the mass of 
the fuel participating in the reaction. The mass of moisture is 
then m H2 o = uiH 2 o m F and the mass of dry fuel is 
m DF = (1 - )tn F . Dry fuel can be further divided into ash and 

dry, ash-free fuel; the mass of ash is m Ash = h m DF = u;” h (l - 
u, H 2 o) m F and the mass of dry, ash-free fuel is 
m D AF = (1 - tUAsh) m DF = (1 - t^Ash)! 1 “ UJ H 2 o) m F- Thus, the ratio 
of the mass of moisture to the mass of dry, ash-free fuel is 


the actual value of the last term in the denominator of Equation 
(27) will be approximately 0.5 ■ 0.99 • 2.4 • (-0.1) = -0.119, which 
means that the actual value of the denominator will be 
1 + 0.045 - (-0.119) = 1.164. Thus, the error will be +5%. A 
similar calculation could be carried out for the carbon dioxide 
generation method. 

7.3. Estimation of the heat release rate when the effect of 
moisture is included 


m H 2 Q _ W H 2 Q _ 1 / gl % 

m DAF l-ui^o 1 - U)" h V 

However, the ratio of the mass of moisture to the mass of 
dry, ash-free fuel can also be computed directly from Equation 
(75). Thus, 


m H 2 0 _ 7 Mh 2 o 

m DAF y Me + 2mM H + nMo 
and we obtain 


(82) 


7 _ M c + 2mM H + nM 0 1 

y“ Mh^o 1 - >"h 2 o 1 - u4*sh 

It can be seen that the effect of ash content on y/y is very 
weak unless exceptionally high values of ash content are 
encountered. In most cases, it can safely be assumed that 
= °- 

Fig. 3 shows the value of m eff + n eff as a function of mois¬ 
ture content for the hypothetical wood described in Section 
7.1 (and for which it is assumed that = 0). It can be seen 
that the effect of moisture content is fairly strong especially 
after the moisture content exceeds approximately 30%. 
However, the effect of this point on the overall accuracy of the 
oxygen consumption method and the carbon dioxide gener¬ 
ation method is fairly small as will be illustrated in the 
following. 

Assume that we set m + n = 13 in the computations (since 
this is the average value for dry wood). Also assume that 
Axo 2 = 0.1. As was discussed in Section 6.4, the value of the 
denominator of Equation (27) in the computations will then be 
1.11. Now, assume that the moisture content of the firewood is 
30%; thus, the value of m e[f + n eff willbe approximately 2.4. Thus, 


The equations for computing the heat release rate were 
derived in Section 6, but those equations were only valid for 
the combustion of dry fuel. Now the equations must be 
modified in two ways. 

Firstly, we want to account for the fact that the moisture in 
the fuel decreases the heating value. Thus, Q(( F in Equation 
(27) is replaced with Qf and Q£' F in Equation (37) is replaced 
with Qf'. Procedures on how to compute the values of Qf and 
Q^" were presented in Section 7.1. 

Secondly, we want to account for the fact that the moisture 
in the fuel changes the composition of the flue gas. Thus, m 
and n in Equations (27) and (37) are replaced with m eff and n eff . 
Procedures on how to compute the values of m eff and m eff were 
presented in Section 7.2. 

Thus, the basic equation of the oxygen consumption 
method becomes 


Q = —NwFGfeQr • 


1 + y (lUeff - tl e ff) 


Ax, 


O 2 


1 + 2 ( m eff ” n eff) (l-*o 2 


-fe(m eff + n eff ) Axo 2 
(84) 


and the basic equation of the carbon dioxide generation 
method becomes 


A y 

Q = NwFGfeQ^"-j-^- (85) 

1 + 2 ( m eff - tleff) *C0 2 + 2 fe ( m eff + tl e ff) AX C o 2 

Theoretically, these equations are valid for fuels containing 
any amount of moisture. As the moisture content approaches 
zero, Equations (84) and (85) reduce to the Equations (27) and 
(37) which were developed for dry fuel. 



Fig. 3 - The effect of fuel moisture on the value of 
m eff + n ef[ for the hypothetical wood described in Section 
7.1. [ ] means that a quantity is dimensionless (e.g. a mole 
fraction). 


8. Uncertainty analysis 

Several potential sources of uncertainty have been identified 
above. In this section, a brief summary is presented. 

8.1. Uncertainty in parameter values 

The uncertainty level associated with the values of the 
parameters Qf, Qf', m eff and n eff is determined by the uncer¬ 
tainty level in fuel characterization. Lowest uncertainty is 
expected if the heating value, moisture content, and 
elementary analysis of the fuel are all known. This case is 
designated as Cl (Complete Information). In this case, accurate 
values for the parameters can be employed, and the only 
remaining uncertainty is caused by the fact that the chemical 
composition of volatile matter (which mainly participates in 
combustion reactions in the early stages of burning) differs 

















18 


BIOMASS AND BIOENERGY 6 I (2OI4) I-2 4 


from the chemical composition of char (which mainly 
participates in combustion reactions in the later stages of 
burning). In the case of wood firing, the magnitude of this 
uncertainty is expected to be of the order of ±5%. It should be 
noted that this uncertainty will remain even if the experi¬ 
mental data include both oxygen concentration and carbon 
dioxide concentration. 

Somewhat higher uncertainty is expected if the heating 
value and moisture content are known but the elementary 
analysis is not known. This case is designated as TI (Typical 
Information). In this case, the values of Q£, Q£", m eff and n e ff are 
partly based on typical values of m and n and partly on actual 
values of heating value and moisture content. Since the values 
of m and n are confined to a fairly narrow range regardless of 
wood species, the uncertainty is not expected to increase 
significantly when compared with the situation where the 
actual values of m and n can be solved from the elementary 
analysis of the fuel. 

If no fuel analyses are available, the uncertainty will be 
considerable. This case is designated as NI (No Information). 
In this case, the uncertainty level in the values of Q£ F and Q£' F 
is of the order of ±10%. In addition to that, any moisture in the 
fuel will have an effect on the values of Q£', Qj,", m eff and n eff , as 
was discussed in Section 7, and the magnitude of this effect 
cannot be estimated unless we have at least some idea on the 
moisture content of the fuel. As an example, consider the case 
where the computations are carried out by assuming that the 
fuel is completely dry while the actual value of the moisture in 
the fuel is 30% It was shown in Section 7 that in this case, the 
effect of the moisture on the heating value of the fuel will 
introduce an error of ±5% in the computed values of the heat 
release rate and that the effect of moisture on the composition 
of the flue gas will cause an additional error of ±5% in the 
computed values of the heat release rate. These two errors are 
additive, so the total error will be ±10%. 

For practical work, we propose the following solution: take 
the typical values Q£ F = 425 kj/mol, Q£' F = 446 kj/mol, m = 0.7 
and n = 0.6 and assume that the moisture content of the fuel is 
15%. For this set of parameter values, we obtain 
Qj? = 415 kj/mol, Qj?' = 436 kj/mol, m e ff = 0.925 and n eff = 0.825. 
Using these values in the computations means that the total 
uncertainty caused by the fact that we do not know the 
moisture content of the fuel will be ±5% as long as the mois¬ 
ture content of the fuel does not exceed 30%. Of course, one 
needs to be careful about the moisture content. The moisture 
content of freshly cut wood easily exceeds 30%, but the wood 
used as firewood is normally dried in air for some period of 
time to achieve a moisture content below 30%. Thus, the 
solution proposed here appears to be reasonable in most but 
perhaps not in all cases. 

8.2. Uncertainty in experimental measurements 

During the development of the oxygen consumption method 
for use in fire safety studies, it was very quickly realized that 
attention must be paid to the uncertainties related to the 
experimental measurements of flue gas flow rate and 
composition. In particular, it was noted that the oxygen 
analyzer must be extremely accurate [6]. However, it will not 
be sufficient just to have a good oxygen analyzer. As the 


method is essentially based on the measurement of the 
difference in the oxygen mole fraction in flue gas and 
combustion air, the operating conditions in the experimental 
setup must be chosen in such a way that the difference is 
substantial enough to achieve good signal-to-noise ratio. 
Thus, the volume flow rate in the exhaust duct must not be 
excessively high. Fitting these two requirements together can 
be rather challenging in some experimental situations. 

Yeager [27] has presented a comprehensive analysis of one 
particular experimental setup, and concluded that an uncer¬ 
tainty level of ±6% could be achieved over a wide range of heat 
release rates provided that the operating conditions were 
chosen appropriately. This level of uncertainty also includes 
the uncertainty in the value of the quantity of heat released 
per quantity of oxygen consumed in the reactions, which 
Yeager stated to be ±5%. Assuming (as Yeager did) that all 
uncertainties were independent of each other, it can be back- 
calculated that his estimate of the combined uncertainties in 
the experimental measurements was approximately ±3%. 

International standards for fire testing usually specify that 
when test methods utilizing the oxygen consumption method 
are calibrated using pure gaseous and liquid fuels with known 
compositions and heating values, the measured values of heat 
release rate must be within ±5% of actual values [14,13], This 
requirement agrees well with the results achieved by Yeager. 

8.3. Total uncertainty 

If the heating value, moisture content and chemical compo¬ 
sition of the fuel are known, then the total uncertainty in heat 
release rate calculations will be approximately ±7%. This 
uncertainty is caused by our not knowing whether the vola¬ 
tiles or the char or both are burning at any given instant of 
time and by the uncertainties related to the experimental 
measurements. 

If the heating value, moisture content and chemical 
composition of the fuel are not known, then the total uncer¬ 
tainty in heat release rate calculations will increase to 
approximately ±13%. 

8.4. Summary 

Based on discussions presented in Section 6, we conclude that 
both the number of moles of oxygen consumed in the 
combustion reaction of wood and the number of moles of 
carbon dioxide generated in the combustion reaction of wood 
can be computed from experimental data obtained using an 
“open system” type setup as described in Section 5. Thus, both 
the oxygen consumption method and the carbon dioxide 
generation method can be used to compute the heat release 
rate during experiments on wood firing in stoves. Normally, 
the equations to be used would be Equation (84) for the oxygen 
consumption method and Equation (85) the carbon dioxide 
generation method. 

To utilize the oxygen consumption method, an oxygen 
analyzer is needed but a carbon dioxide analyzer is not 
necessary; to utilize the carbon dioxide generation method, 
a carbon dioxide analyzer is needed but an oxygen analyzer is 
not necessary. As the prices of gas analyzers have become 
more affordable, it is nowadays a common practice to use 
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both an oxygen analyzer and a carbon dioxide analyzer for 
flue gas analysis. For heat release rate analyses, this provides 
an additional way of checking the results. 

It is fully realistic to assume that the uncertainty level in 
the values of the heat release rate can be pressed down to ±7% 
if the experimental setup is appropriately implemented and 
the properties of the fuel are comprehensively analyzed. 


9. Application to an actual stove experiment 

9.1. Test room 

A test room for studying the performance of heat storing 
stoves was built inside a laboratory hall at VTT. The objective 
was to provide means for accurate determination of the heat 
output from the stove as a function of time. This was quite 
challenging since the heat output varies considerably with 
time and the longest tests take several days. 

The solution adopted at VTT was based on the concept of 
keeping the test room at constant temperature throughout the 
whole test. The test room was lightly structured, well insu¬ 
lated and equipped with a controllable cooling system 
designed to maintain a constant temperature inside the test 
room. The air temperatures inside and outside the test room 
were equal; thus, the uncertainty associated with the tran¬ 
sient behavior of the walls could be eliminated. During the 
tests, the heat flow extracted by the cooling system was 
measured continuously and discharged outside the labora¬ 
tory. Since there was no heat accumulation into the test room 
walls, the heat flow extracted by the cooling system was equal 
to the heat output of the stove. 

It has to be noted, however, that the heat output from the 
stove to the test room and from the test room to the 
surroundings are of no importance to the analyses presented 
in this paper. We are only considering the heat release rate 
from the burning firewood inside the stove, and this 
measurement does not at all depend upon the temperatures 
of the stove or the test room. This is also the reason why we 
have chosen to use the term 'oxygen consumption method’ 
and not the commonly used term ‘oxygen consumption 
calorimetry’. Of course, the way the energy released from the 
burning firewood is stored in the stove structure and trans¬ 
ferred from the stove to the test room is quite important from 
the point of view of thermal comfort and energy efficiency of 
the whole system, but we do not get into those aspects here. 

Combustion air was supplied from outside the test room 
directly to the stove through a duct which was connected to 
a windbox attached to the stove. The stove was connected to 
a chimney, the top of which was located underneath an 
exhaust hood. An induced draft fan was employed in the hood 
to enhance the draft in the chimney; such an arrangement 
was needed in order to emulate realistic situations where the 
chimney is taller and the cold temperature of outside air helps 
to create a stronger natural draft. An additional advantage 
was that using a fan ensured that the flue gas exiting from the 
chimney could not escape elsewhere but was extracted 
through the hood. The induced draft fan was in operation 
throughout the whole experiment; thus there was a small flow 
rate through the chimney already before the test was started. 


When the firing period ended and the chimney damper was 
closed, the flow rate through the chimney decreased to zero. 

During the experiments, the temperature and flow rate of 
combustion air and the temperature, composition and flow 
rate of the flue gas through the chimney were continuously 
measured (see Section 9.3). It should be noted, however, that 
the walls of the test room were not gas-tight and some of the 
combustion air therefore flowed into the room through leaks 
in the walls. As has been discussed previously, the existence 
of such leaks makes it more attractive to determine the heat 
release rate by using the open system arrangement (in which 
the measured data on the flue gas flow rate is used in the 
calculations) rather than by using the closed system 
arrangement (in which the measured data on the combustion 
air flow rate is used in the calculations). This is also what we 
did, as will be described in Section 10. 

9.2. Stone 

The stove was of the type Saara made by Kerman Savi in 
Finland. The height of the stove was 2100 mm and the mass 
was 1350 kg. The stove is shown in Fig. 4. 

9.3. Instrumentation 

The flow rate of combustion air was measured using an Alnor 
Thermo-Anemometer Type GGA-26 Serial No 850-088 and the 
relative humidity of combustion air was measured using 
a Rotronic Hygroskop manufactured by Fattore Vitale. 

The flow rate of flue gas in the chimney was measured 
using a calibrated orifice plate equipped with measurements 
of pressure drop, pressure and temperature. The mole 



Fig. 4 - Stove installed inside the test room. Ductwork and 
windbox for introducing the combustion air and 
thermocouple wires for stove temperature measurements 
are seen in the figure. 
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fractions of oxygen, carbon dioxide and carbon monoxide in 
dry flue gas were measured with IPA gas analyzers manufac¬ 
tured by PPM Systems. The ranges were 0—20% for C0 2 and 
0-25% for 0 2 . The precision of the analyzers was ±2% of full 
scale. Flue gas sampling for analysis (temperature and 
concentrations) was arranged in accordance with the 
following standards: prEN15250, prEN13240 (fireplaces), EN304 
and EN303-5 (boilers). 

The measurement of the flue gas flow rate provides infor¬ 
mation essentially in real time, whereas the measurement of 
flue gas composition is inevitably accompanied by a consid¬ 
erable time lag. The time lag is caused both by the delay time 
(the time which the gas sample needs for traveling from the 
sampling probe located in the chimney to the gas analyzers 
located in an instrumentation rack at some distance from the 
chimney) and by the response time (the time which the 
analyzers need to react to a change in the composition of the 
sample gas). In fire testing, considerable attention is being 
paid to the delay and response times and it is normal practice 
to apply time shifting of the data to achieve synchronization 
(see, e.g., [14]). This is very appropriate considering how 
rapidly natural fires develop (see, e.g., the plots of heat release 
rate vs. time in the cases compiled by Babrauskas [2]). In our 
case, events developed much more slowly and no attempt was 
made to compensate for the time lags of the gas analyzers. 

It should be noted that of the measurements listed above, 
only the relative humidity of combustion air, the flow rate of 
flue gas in the chimney, and the mole fractions of oxygen and 
carbon dioxide in the flue gas are needed for the determina¬ 
tion of the heat release rate by the open system arrangement 
employed in this study. 

Considerable number of other measurements were also 
carried out during the test; e.g., the internal temperature 
distribution of the stove was recorded by using thermocouples 
embedded in the structure of the stove during erection. Those 
measurements were not needed for the heat release rate 
analyses presented in this paper but were very useful for the 
purposes of development and optimization of the stove. 

9.4. Fuel 

The fuel consisted of 50 mm x 50 mm x 330 mm sticks made 
of Finnish softwood. The properties of the fuel are presented 
in Table 12. The total amount of fuel used in the test was 
8.97 kg corresponding to an energy content of 143.3 MJ. 


For this fuel, the chemical composition corresponds to 
m = 0.744 and n = 0.605 (see Section 4.1). Thus, m + n = 1.349 
and m - n = 0.139, which values agree well with the values 
typically found for different wood species. 

The heating value of the dry fuel is 18.9 MJ/kg, which can be 
inserted to Equation (9) to yield Q£ F = 411 kj/mol. This is the 
amount of heat released during combustion of dry wood 
divided by the number of oxygen moles consumed in the 
combustion reactions. Similarly, Equation (11) yields 
Q£'p = 440kJ/mol. This is the amount of heat released during 
combustion of dry wood divided by the number of carbon 
dioxide moles generated in the combustion reactions. The 
values thus obtained are slightly lower than the generic values 
425 kj/mol and 446 kj/mol but well within the ±10% uncer¬ 
tainty range postulated by us. 

Further, Equations (70), (72) and (74) can be used to evaluate 
the effect of fuel moisture. Thus, we obtain Q(' = 403 kj/mol and 
Q(7 — 431 kj/mol. These results confirm the conclusion 
expressed in Section 7.1 that in the context of using the oxygen 
consumption method and the carbon dioxide method, the effect 
of fuel moisture on the heating values is indeed very small 
except at fairly high values of moisture content; in our case, 
a moisture content of 13.7% by mass only decreased the heating 
values by 2% when compared with the values for dry wood. 

Evaporation of fuel moisture during the combustion not 
only represents an energy loss (which is reflected in the values 
of Q(' and Q("), but also has an effect on the composition of the 
flue gas. To take this effect into account, we replace m and n 
with m e ff and n e ff as suggested in Section 7.2; for our fuel, we 
have meff = 0.949 and n e ff = 0.810. 

As for parameter £, Equation (42) predicts a value of 0.948, 
which can be used to check the consistency of the gas 
analyzer outputs. The value of £ is independent of the mois¬ 
ture content of the fuel. 

9.5. Experimental procedure 

The test described here was carried out in August 2002. The 
firing of the stove was started at 10:20 a.m. on 26 August and 
was completed at 12:05 p.m. on the same day. The monitoring 
of the temperatures within the stove structure and the 
measurement of the heat output from the stove to the test 
room continued till 28 August, but for the purposes of the 
present paper, we are only interested in the firing period 
which lasted for the first 105 min of the test. 

The events during the test are listed in Table 13. 


Table 13 - Events during the test. 

Time 

Elapsed time 

Event 

h:min 

min 


10:20 

0 

Insertion of the first batch of fuel (3.50 kg) 
and ignition 

10:55 

35 

Insertion of the second batch of fuel 
(2.98 kg) 

11:27 

67 

Insertion of the third batch of fuel (2.49 kg) 

11:54 

94 

Stirring of embers on the grate, draft 
decreased 

12:05 

105 

Combustion air damper and flue gas 
damper closed 


Table 12 - Properties of the fuel used in the test. 

Property 


Value 

Method of 
determination 

Moisture content 

As fired 

13.70% 

DIN 51718 

Lower heating value 

As fired 

15.98 MJ/kg 

DIN 51900 

Lower heating value 

Dry basis 

18.90 MJ/kg 

DIN 51900 

Ash content 

Dry basis 

0.2% 

DIN 51719 

Carbon content 

Dry basis 

51.6% 


Hydrogen content 

Dry basis 

6.4% 


Nitrogen content 

Dry basis 

0.2% 


Oxygen content 

Dry basis 

41.6% 

(by balance) 
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Time [h:min] 

Fig. 5 - Flue gas flow rate during the test. 


10. Results 

10.1. Flue gas flow rate and composition 

Flue gas flow rate in the chimney during the test is shown in 
Fig. 5. The insertion of the second and third batches of fire¬ 
wood had little effect on the flue gas flow rate, whereas the 
decreasing of the draft at 11:54 and the closing of the chimney 
damper at 12:05 can be clearly seen. 

Oxygen mole fraction in dry flue gas during the test is 
shown in Fig. 6, carbon dioxide mole fraction in Fig. 7 and 
carbon monoxide mole fraction in Fig. 8. The mole fractions of 
oxygen and carbon dioxide are almost like mirror images of 
each other, as was expected. Firewood insertions and 
decreasing of the draft are clearly seen in the data. 

Carbon monoxide mole fraction is quite low during the 
whole test, typically only a few percent of the mole fraction of 
carbon dioxide. Thus, CO may be safely ignored in heat release 
rate analysis [28] (in such a situation, the effect of CO on heat 
release rate analysis is well below 1%). The total yield of 
carbon monoxide during the whole test is approximately 29 g 
CO per 1 kg dry mass of the fuel. 

The value of parameter £ is shown in Fig. 9. For most of 
the time, the value of £ varied between 0.96 and 0.98, which is 
in good agreement with the expectations based on the 



10:00 11:00 12:00 13:00 

Time [h:min] 

Fig. 6 - Oxygen mole fraction in dry flue gas during the 
test. 



Time [h:min] 


Fig. 7 - Carbon dioxide mole fraction in dry flue gas during 
the test. 


composition of the fuel. Randomly fluctuating values were 
observed at the very beginning of the test, when the gas 
compositions had not yet really changed, and at the end of 
the test, when the gas compositions had already returned 
back to the normal values for atmospheric air. At those 
situations even the smallest errors in experimental 
measurements will have a large effect on the value of £ and 
we therefore should not expect reliable results at these 
conditions. Thus, the gas analyzers also appear to have been 
working properly. 

10.2. Heat release rate by the oxygen consumption 
method 

In the following, we will present results computed in three 
different ways. The baseline case is designated as OX-NI, 
where NI stands for “no information”. This is a hypothetical 
case, where we assume that nothing at all is known about the 
fuel, and the recommended values Q£ F = 425 kj/mol, m = 0.7 
and n = 0.6 are used in the calculations. Further, the fuel is 
assumed to be completely dry. These assumptions mean that 
the heat release rate will be overestimated (both because the 
moisture content represents an energy loss and because the 
evaporation of the moisture adds ballast to the flue gas flow 
rate). It is clear that with such assumptions the accuracy 



Time [h:min] 

Fig. 8 - Carbon monoxide mole fraction in dry flue gas 
during the test. 
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Fig. 9 - The value of parameter £ during the test. 


cannot be expected to be particularly good, but the baseline 
case is nevertheless interesting as an indication of what can 
be achieved in the worst possible conditions. 

The second case is designated as OX-TI, where TI stands 
for “typical information”. Here we assume that the heating 
value and the moisture content of the fuel are known, but the 
elementary analysis is not. This would probably be quite 
typical in the situations where the oxygen consumption 
method would be applied. Thus, we will use the actual values 
of the heating value and moisture content and the recom¬ 
mended values m = 0.7 and n = 0.6 to obtain values for Qf, m eff 
and n eff . The accuracy is expected to be much better than for 
the baseline case since the variations in the values of the 
parameters m and n are fairly small for most wood species of 
interest. 

The third case is designated as OX-CI, where Cl stands for 
“complete information”. Here we assume that we know not 
only the heating value and the moisture content of the fuel, 
but also the elementary analysis. Of course, such a situation 
will provide the means to achieve the best possible accuracy. 
In the present study we indeed do have the complete 
information of the fuel, but normally this cannot be 
expected. 

The results for all three cases are shown in Fig. 10. It can be 
observed that somewhat higher values for the heat release 
rate are obtained for the case NI than for the other two 
cases. For the cases TI and Cl, the results are almost undis- 
tinguishable. All this is quite in agreement with what was 
predicted earlier. 

10.3. Heat release rate by the carbon dioxide generation 
method 

Similarly to the oxygen consumption methods, we will 
present results computed in three different ways: NI, TI and 
CL For the NI case, we set Q£' F = 446 kj/mol, m = 0.7 and n = 0.6 
and assume that the fuel is completely dry. 

The results for all three cases are shown in Fig. 11. A 
pronounced resemblance to Fig. 10 can be observed. Again, 
case NI produces somewhat higher values for the heat release 
rate than the other two cases. The results for cases TI and Cl 
are almost undistinguishable. All this is quite in agreement 
with what was predicted earlier. 



Time [h:min] 


Fig. 10 - Heat release rate during the test as computed 
using the oxygen consumption method. The cases TI 
(typical information) and Cl (complete information) are 
almost undistinguishable. The case NI (no information) 
predicts somewhat higher values for the heat release rate 
as was anticipated in the text. 


10.4. Comparison of heat release rate analyses 

The results shown in Figs 10 and 11 can be integrated with 
time to obtain the total amount of heat released during the 
test. The results are compiled in Table 14. 

It can be seen that both the oxygen consumption method 
and the carbon dioxide generation method seem to be 
working very well. For the TI and Cl cases, the difference 
between the calculated results and the actual value is only 
a few percent; better agreement can hardly be expected 
considering what experimental uncertainties must be dealt 
with. For the NI case, the discrepancy is lower than 10% which 
is still quite satisfactory. 

Finally, it may be mentioned that the total amount of 
carbon dioxide generated during the test is 14.65 kg when 
computed from the composition and total amount of wood 
inserted into the stove during the test and 15.12 kg when 



Time [h:min] 


Fig. 11 — Heat release rate during the test as computed 
using the carbon dioxide generation method. The cases TI 
(typical information) and Cl (complete information) are 
almost undistinguishable. The case NI (no information) 
predicts somewhat higher values for the heat release rate 
as was anticipated in the text. 
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Table 14 - Summary of calculated values for the total heat 
release. The calculated values are compared with 
143.3 MJ, which is the total energy content of the fuel 
inserted into the stove during the test. 


Case 

Total heat 
release 

Difference 

Case 

Total heat 
release 

Difference 

OX-NI 

152.8 MJ 

+6.6% 

CA-NI 

155.2 MJ 

+8.3% 

OX-TI 

145.6 MJ 

+1.6% 

CA-TI 

146.8 MJ 

+2.4% 

OX-CI 

144.7 MJ 

+1.0% 

CA-CI 

148.5 MJ 

+3.6% 


computed from the flue gas flow rate and mole fraction of 
carbon dioxide in the flue gas. Thus, the combined error of 
experimental measurements is just slightly above 3%, which 
also represents fairly good quality of experimental work. 


11. Conclusions 

Both the oxygen consumption method and the carbon dioxide 
generation method proved to be excellent methods for the 
determination of the heat release rate during the combustion 
of firewood in stoves. The accuracy of the methods appears to 
be very good, particularly if detailed information on the 
composition and heating value of the fuel are available. The 
computations are somewhat lengthy, but can be easily 
implemented on a computer spreadsheet program. 

The methods are independent from each other in the sense 
that the oxygen consumption method does not require 
information on the carbon dioxide content in the flue gas and 
the carbon dioxide method does not require information on 
the oxygen content in the flue gas. Thus, the experimental 
setup can be kept fairly simple. Since it is nowadays 
a common practice to use both an oxygen analyzer and 
a carbon dioxide analyzer for monitoring the combustion 
process, the oxygen consumption method and the carbon 
dioxide generation method can be both applied to the same 
experiment to provide a mutual self-check. Of course, it has to 
be kept in mind that both methods use the same information 
of the flue gas flow rate and errors in flue gas flow rate 
measurement may therefore pass this self-check undetected. 
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